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ABSTRACT 


During the summer and fall of 1970, a line of 
nine magnetic observatories was set up along a corrected 
geomagnetic meridian through western Canada. This thesis 
contains an analysis of polar magnetic substorms recorded 
at these stations and at a large number of standard obser- 
vatories around the world. The research was directed 
toward the study of the dynamic development of the substorm 
current system, and the accurate modelling of the current 
system at various instants throughout the substorm. Data 
from all-sky cameras on the station line facilitated a 
study between substorm current system fluctuations and the 
variations in auroral arc structure. 

It was found that the substorm current system 
develops through the growth of the westward electrojet in 
the midnight sector of the auroral oval. The substorm 
further develops through quasi-periodic intensifications 
of the northern border of the electrojet in conjunction 
with the generation of irregular loops and surges in the 
auroral arc structure. An intensification of the eastward 
electrojet in the evening sector accompanies the substorm 
onset, although the eastward electrojet may be active with- 


out accompanying substorm activity. 
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Model studies of the electrojet yielded several 
plausible three-dimensional current system configurations 
which could explain the substorm perturbation patterns. 
Distortions of the northern border of the electrojet were 
generally found to be represented by confined north-to- 
south current flow which was capable of explaining the 
large positive D-component perturbations in the region of 
the surges or loops. It was found to be extremely impor- 
tant to include the effects of induction in the modelling 
of the current systems, and excellent agreement between 
the perturbation patterns by the model current systems and 
the observed perturbation patterns resulted when induction 
was taken into account. 

Several different methods of portraying the 
dynamic development of the auroral electrojets are employed, 
and it is found that no one presentation is completely 


adequate to transfer all the information involved. 
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CHAPTER I INTRODUCTION 


1.1 The magnetospheric substorm 


It has long been realized that the energy which is 
dissipated in high latitude auroral and electromagnetic 
disturbances is derived from the interplanetary medium 
through the interaction of the solar wind with the outer 
regions of the Earth's magnetic field configuration. In 
this interaction, the Earth's dipole magnetic field is 
highly distorted, assuming a configuration as shown in 
Fig. 1.1. The magnetic field configuration, known as the 
magnetosphere, features a compression of the dipole lines 
on the dayside and a distention of the magnetic lines of 
force on the nightside forming the magnetotail. The 
magnetic field pattern shown in Fig. 1.1 is representative 
of steady state conditions, and the action of the solar 
wind on this configuration results in a flow of energy 
into the inner magnetosphere. While there are conflicting 
theories regarding the exact fashion in which energy may 
penetrate into the inner magnetosphere (eg. Axford and Hines, 
1961; Dungey, 1961; Frank, 1971), it has become apparent 
that the rate of energy transfer is highly dependent on the 


parameters of the interplanetary plasma, and in particular 


et dotdw yprons ont tsd¢ bostisor aed pnp t ean bt 
ottenpsmoytosls bas {svo1us sbuttiel dptd nt bsteqteath- 
murbem yrstonsfarstai ont mort bevinsb 2t seonsdvut ath 
netuo elt dtiw batw vslo2e sit to wotfosysdat ont igor 

nl .softstsywettnos blott asttanpsm 2'd3153 eat To znoteot 

et blot? attenpsm sfogib 2'di1s3 edt _nottosrstat ate 

nt nwode 26 nofssyupitroo 6 patmuees ebotrotath vido td 

ait 26 nwondt ,notdesupttncs blatt sitengem sAT tat wae 
2anrl sfoqrb sit to moteesyqmos 5 es yussat 9194920) m 

to 2ontl oitenpem aht-to nokinotetb 5 bns ohfeysb ond no 

ont .ffstotonpsm sit pnimvot sbtetdptn sit no so40F 
avitsinsesiqey et [.f .pri nt nwode nvettsq biat? ot tenpem 
wsloe oft to nortos sit bns5 .2norstbnos ststz2 ybsst2 to 
verens to wolt 6 nt 2etivest notte vupttnos 2d? ao batw 

| enttotltnos a6 sont af KAW -sysigzotonpem tenant ady osnt 
Vem verene fofdw nt notdest tosxe sit ontbyepey estrosdd 
ezenth bas byotxA .p9) evedqzotenpsm tenai st osnt sisiseneg 






ec ete Sr nse 0 om oneal 
. es ail : a : 
. Pe | 


_ nega enooed a oF TG eafier7 :faer , haagaes. 2 _ 















7 





*aueitd LeUuoLpLuow 


qUBLUpLW-UOOU ILUaYydsozZaUBew ueLOS ayy UL BuUaydsozyaubew ay} 40 uoLZeunBLyuod abeuaray |{°| “bL4 









































ASNVdOLANSVA 











ASNVdVWSV Id 





sdetaibim-noon sirtsigeotenpsm ysloe odd nf syanqeotenesm sit to notte ywertaes 
.onelq Ts 





the direction of the component of the interplanetary magnetic 
field normal to the ecliptic plane. Regardless of the means 
by which the energy enters the magnetosphere, it appears 
to be stored in the magnetotail through the distention of 
LG Leh T Tetra trite ss 

While energy may be dissipated slowly and 
continually within the magnetosphere through the generation 
of quiet auroral arcs, ionospheric current systems etc., 
the most outstanding magnetospheric activity is caused by 
sudden impulsive releases of energy. It is generally felt 
that the energy is released through explosive magnetic field 
lines merging in the magnetotail (Atkinson, 1966; Dessler, 
1968) and during such periods of release of energy, many 
diverse electromagnetic phenomena are generated. The sum 
of all these impulsive phenomena has been termed the 
magnetospheric substorm, and the character of this pertur- 
bation has been discussed in great detail by Akasofu (1968). 
He has subdivided the magnetospheric substorm into substorms 
involving different portions of the electromagnetic spectrum, 
viz. the auroral substorm, the polar magnetic substorm, the 
ionospheric substorm, the X-ray substorm, the proton aurora 
substorm, the VLF emission substorm, and the micropulsation 


substorm . Throughout this thesis we will be primarily 
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interested in the auroral and polar magnetic substorms. For 
a comprehensive review of all the various substorm phenomena, 


the reader is referred to Akasofu (1968). 


1.2 The auroral and polar magnetic substorms 


As we have stated in the previous section, the 
auroral and magnetic signatures of the magnetospheric substorm 
are termed the auroral substorm and polar magnetic substorm 
respectively. At present the development of these two 
substorms is best described in the framework developed for 
the auroral substorm by Akasofu (1964) with some modifica- 
tions involving more recent observations. 

The development of the auroral substorm has been 
subdivided into the expansive phase and the recovery phase, 
and each of the phases has been divided into three stages 
as follows (see Fig. 1.2 and 1.3). 

(i) Expansive Phase Stage I: The auroral arcs 
near midnight brighten suddenly, at which time the magnetic 
field changes sharply in the vicinity of the arc and a Pi 2 
micropulsation burst is initiated. 

Stage II: The auroral arcs 
become distorted and move poleward, during which time the 


magnetic field perturbation tends to intensify. 
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Fig. 1.2 Examples of a polar magnetic substorm recorded at a 
high-latitude station (Fort Churchill) and a mid- 
latitude station (Victoria). Note the burst of 
Pi 2 micropulsation activity associated with the 
substorm onset. A typical polar magnetic substorm 
in the high-latitude regions has an amplitude on 
the order of a couple hundred gammas (one gamma 
equals 107° Gauss). 
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Schematic diagram to illustrate the development of 
the auroral substorm. The center of the concentric 
circles in each stage is the north geomagnetic 
pole, and the Sun is toward the top of the diagram 
(after Akasofu, 1964). 
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Stage III: The westward 
travelling surge and eastward drifting patches are formed, 
at which time the magnetic field perturbation intensifies 
eon. V.: 

(ii) Recovery Phase Stage I: The arcs reach their 
northernmost position and start to recede southward, while 
Surges and loops distort the northern border. During this 
time, the magnetic field perturbations remain large and 
irregular. 

Stage II: The arcs continue 
to drift equatorward and the drifting loops change into 
irregular bands. The magnetic perturbations start to die 
down at this time. 

Stage III: The arcs slowly 
complete their equatorward drift, and the magnetic activity 
ceases. 

It has recently become apparent that there is a 
period of distinctive activity in the magnetosphere prior to 
the onset of the expansive phase of the substorm. This pre- 
substorm activity has been termed the growth phase (McPherron, 
1970), and is thought to represent the period of time during 
which there is an enhanced flow of energy from the solar wind 


into the magnetosphere. At the present time there is consider- 
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able controversy regarding the auroral and ground-based 
Magnetic signatures of the growth phase. While McPherron 
(1970) claims that a Slowly deepening H-component trace at 
auroral zone observatories is typical of the growth phase, 
this particular perturbation is also the Signature of a 
substorm expansive phase seen some distance from the intensi- 
fied portion of the auroral electrojet. In addition Feldstein 
(1972) and Mozer (1971) have Suggested that equatorward 
drifting arcs signify the development of the growth phase, 
while Akasofu and Snyder (1972) and Subbarao and Rostoker 
(1972) claim that this is not the case. In this thesis we 
will comment in more detail on the magnetic signatures of 


the various phases of substorm activity. 


1.3 Current systems believed to be responsible for polar 
a A SS ee? phe teS>R0 RF) O1L0 phOF apa Lak 


magnetic substorms 


For some time it has been known that current systems 
in the Earth's upper atmosphere are responsible for the large 
magnetic variations observed during auroral displays (see 
Chapman, 1967). The first definitive study of such current 
systems was done by Birkeland (1908, 1913), who concluded that 
the current system associated with, what he termed “Otay 


elementary storm', was three-dimensional in nature with current 
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flowing down along magnetic fields to the ionosphere, 
westward through the ionosphere and back up the field lines 
into space. He also introduced the equivalent current 
concept for the purpose of data presentation. Chapman (1935) 
disagreed with Birkeland's model and suggested that all 
Currents responsible for the observed magnetic perturbations 
were confined to the ionosphere. Vestine and Chapman (1938) 
modelled the auroral electrojet with this idea in mind. 

The next important study of polar substorms was 
made by Fukushima (1953) in which he effectively used the 
equivalent current system to analyze the development of 
geomagnetic bays. In his model studies he used distributions 
of electric doublets confined to the ionosphere. Heppner 
(1954) discussed the possible ionospheric current systems 
which could explain the diurnal patterns of geomagnetic bays 
observed at College, Alaska. He concluded that a two jet 
System, an eastward electrojet in the evening sector and a 
westward electrojet in the morning sector, could explain these 
observations. 

Akasofu et al. (1965) introduced a new model current 
System that included a strong westward electrojet flowing 
along the auroral oval with return currents flowing across 


the polar cap and also at lower latitudes below the auroral 
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oval. They believed this return current to be responsible 
for the positive H bays in the evening sector as observed 
by Heppner. 

About this time the concept of field-aligned current 
as originally proposed by Birkeland was again considered 
(Bostrém, 1964). This followed from the realization that 
certain regions of the magnetosphere were not subject to 
perfect frozen field conditions as proposed by Alfvén and 
Falthammar (1963) (also see Block, 1967). More recent studies 
by Armstrong and Zmuda (1970), Cloutier et al. (1970), 

Vondrak et al. (1971), and Park and Cloutier (1971) have 
confirmed the existence of Birkeland (field-aligned) currents 
(see Cloutier (1971) for a recent review of such currents). 

Model studies of three-dimensional current systems 
which include field-aligned current linkage of the westward 
electrojet to the outer magnetosphere, have been used recently 
to model substorm perturbations by Akasofu and Meng (1969) 
and Bonnevier et al. (1970). Despite the ability of these 
three-dimensional current systems to generate positive H 
regimes at high and low latitudes, the positive bays in the 
evening sector suggest the existence of an eastward electro- 
jet (Langel and Cain, 1967; Feldstein and Zaitzev, 1968) and 


from the work of Cummings (1966), Cummings and Dessler (1967), 
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and Kamide and Fukushima (1971), this electrojet would appear 
to be associated with the asymmetric ring current in the 
evening sector (see Rostoker (1972) for a comprehensive 


review of polar magnetic substorms). 


1.4 The problems studied in this thesis 


Previous studies of polar magnetic substorms were 
made using magnetic data from a distribution of magnetic 
observatories which, in most cases, were Sparsely scattered 
over the surface of the Earth (eg. Birkeland, 1908; Vestine 
and Chapman, 1938; Silsbee and Vestine, 1942; Fukushima, 

1953; Akasofu and Meng, 1969). Since the ionospheric portion 
of the current system responsible for the polar magnetic 
substorm is limited in both latitudinal and longitudinal 
extent, a detailed analysis of the dynamic development of 
polar magnetic substorms was not possible. This fact prompted 
us to install a line of magnetometer stations along a common 
geomagnetic meridian through Western Canada. This thesis 
pertains to the acquisition and analysis of data obtained 

from this line of stations during the summer of 1970. 

In the following chapters we will present a detailed 
analysis of several polar magnetic substorms, outlining the 


development of the polar electrojets and the associated auroral 
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activity. We will use three-dimensional equivalent current 
Systems to model the substorm current system. We will show 
that when induction and asymmetric current distributions 
are taken into account, excellent agreement between the 
observed magnetic perturbation pattern and the theoretically 


calculated pattern is obtained. 
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CHAPTER II EVALUATION OF THE MAGNETIC FIELD DUE TO 
THREE-DIMENSTONAL CURRENT SYSTEMS 


Por’ Introd vetion 


In this chapter we present a convenient method of 
evaluating the magnetic field associated with three-dimensional 
current systems using Biot-Savart's law. The advantages we 
find in using this particular method over other methods pre- 
sented in the literature are: (1) the ease with which this 
method can be applied to arbitrary current systems in the 
magnetosphere (the assumption being made that the geometries 
of these current systems can be expressed in polar coordinates 
with the origin located at the center of the Earth), (2) the 
ability to include effects of induction due to a conducting 
Earth (although we assume infinite conductivity, the finite 
conductivity case can also be included using this method), 
and (3) the small amount of computer time required to integrate 
over current systems with large spatial dimensions. 

Other methods for obtaining the magnetic effects 
associated with various three-dimensional model current sys- 
tems have been reported in the literature. These magnetic 
effects may be evaluated either by performing numerical 


computations or by using wire models and a search coil (see 
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Bostrém (1968) for a review of this Subject). Recently Meng 
and Akasofu (1969) revised Kirkpatrick's model system (Kirk- 
patrick, 1952) in order to explain magnetic perturbations 
associated with polar magnetic substorms. Kawasaki and 
Akasofu (1971) computed the magnetic effects associated with 
this revised model. Bonnevier et al. (1970) (also see Bostrom, 
1968) used a distribution of infinitisimal magnetic dipoles 
along dipole field lines to obtain the scalar magnetic poten- 
tial and thus the magnetic field associated with field-aligned 
Currents. We also use this method in conjunction with some 
specific model current systems (see section 2.5 and Appendix 
Al). Kamide and Fukushima (1971) used Biot-Savart's law 

to calculate the magnetic perturbations associated with a 
three-dimensional eastward electrojet. Although we also use 
Biot-Savart's law, our mathematical formulation presented in 
this chapter and Appendix Al is such that the differential 
form of Biot-Savart's law can be represented by a single 
matrix which includes all the necessary coordinate transfor- 
mations. In so doing, a large number of terms either cancel 
or become greatly simplified, thereby allowing more accurate 
computations to be made with a significant decrease in compu- 
ter time. Such a formulation also lends itself to a solution 
of the induction problem for three-dimensional current systems 


(see Appendix A2). 
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2@.2 Evaluation of magnetic fields due to three-dimensional 
ee UE TO Tnree-dimensional 


current systems using Biot-Savart's law 


2.2.1 Arbitrary current loop 


The application of Biot-Savart's law to an 
arbitrary current loop external to the earth's surface 
is presented in Appendix A.1. We showed that for a current 
loop (C) (see Fig. Al.1) the total contribution to the jth 


field component at (1285 94,) is 


where J is the current flowing in the loop, k isa 


proportionality factor dependent upon the system of units 


used and tery is an element of the matrix 
0 ro A3798] “Vj A51 dS] 
dc = 3 rA,3ds, (rA,. + ry A3o)ds5 (rA,, rAoo)ds, 
-rA, ods (r A33 ~ rAj5) ds. -(rA53 * rA35)ds4 
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All variables in dc are defined in Appendix Al. The 


field components given by equation 2.1 are related to the 


total magnetic field at (r 285905) by 


(23:) 


Equation 2.1 is indeed very simple to use for 
evaluating the field due to three-dimensional current 
Systems provided, of course, an analytic expression for 
the current path Cy /exists. This equation also lends 
itself to the problem of determining the field due to 
induced currents in a uniformly conducting earth (see section 


2.3 and Appendix A2, section A2.3). 


2.2.2 East-west (E-W) current system 


A model of a three-dimensional E-W current system 
1s) 41 lustra Ged We barge 2 dicen thts system, field-aligned 
(Birkeland) currents flow from the magnetospheric equatorial 
plane down the field lines to the eastern edge of the 
ionospheric westward electrojet. After flowing westward to 
the western edge of the ionospheric electrojet, they flow 
back up the field lines to the magnetospheric equatorial 


plane. The path of current flow is completed by eastward 
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directed currents in the equatorial plane. In model studies 
this is analogous to a decrease in the ring current intensity 
in this particular sector of the equatorial plane. A model 
system with eastward directed ionospheric currents (east- 
ward electrojet) is obtained merely by reversing the direc- 
tion of current flow. 

Since the configuration of the total current system 
is controlled by the configuration of the magnetic dipole 
field, the physical parameters of the total system can be 
represented by those of the ionospheric segment alone. This 
proves to be very useful since the ionospheric segment pro- 
vides the largest contribution to the magnetic perturbations 
observed by our line of stations. 

In the following discussion we refer to Birkeland 
currents at each end of the westward electrojet as eastern 
and western sheet currents. This does not, however, suggest 
that such current flow is confined to a narrow region in 
longitude as the term 'sheet' implies. In fact such current 
flow is probably distributed along the entire length of the 
westward (or eastward) electrojet. An example of a model 
with this particular current distribution is presented later 
in this chapter. For the present development we assume any 


longitudinal distribution of Birkeland currents can be 
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approximated by sheet currents at both ends of the electrojet. 
We also refer to the ionospheric segment as a sheet current. 
Here, we assume the height-integrated ionospheric segment 

can be approximated by a horizontal sheet current at a height 
h above the surface of the earth. 

To facilitate the following mathematical develop- 
ment we use a series of primary current loops such as those 
shown in Fig. 2.2 to represent an integral part of the total 
current system. The two examples (a) and (b) are primary 
loops with ionospheric current flowing E-W and N-S along 
lines of constant latitude and longitude respectively. In 
general, however, the ionospheric segment and thus the 
complete primary loop can have any orientation and shape 
necessary to coincide with the geometry of the total current 
system. 

The latitudinal and longitudinal extent ('width' 
and 'length' in degrees) of the ionospheric sheet current 
for the model in Fig. 2.1 is shown in Fig. 2.3. The lines of 
x's and dots designate incoming and outgoing Birkeland 
currents respectively. The angles o¢ (central meridian) 
and 86 represent the centers of the current system with 
respect to longitude and to latitude. Both parameters prove 
to be very useful when comparing observations of magnetic 


perturbations with those of model current systems. 
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(B) 


Fig. 2.2 (a) E-W primary current loop associated with the 
model current system illustrated in Fig. 2.1. 
(b) N-S primary current loop associated with the 


model current system presented later in this 
section (see Fig. 2.7). 
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Diagram showing the ionospheric sheet current for 
an E-W current system. Also shown is an arbitrary 
distribution of current across the electrojet (see 
section’ 2.5). 
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Any point P on the ionospheric sheet current is 


P(r, 58) 5o)) where 
i te Bee h (a = radius of earth). 


The ionospheric element of the particular primary loop used 
in the evaluation of the field due to the model in Fig. 2.1 
is also shown in Fig. 2.3. Thus, once we have obtained an 
analytical expression for the field of this particular 
primary loop as a function of Oy » the total field due to 
the current system in Fig. 2.1 can be evaluated by integra- 
ting this expression over the width By to 86 
To allow for arbitrary latitudinal distributions 

of current flow in the system, we define a latitudinal 
current distribution function (6, 59) such that the total 


current J is given by 


Jo.) = 1(0,56,) ry do, (2.5) 


0. (4,) 
8 y (0) 


where we have allowed for a variable width by denoting ors 
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J(,) flowing across any meridian between oe and oy 
is constant and therefore independent of oy » we can 
arbitrarily define both I(6,) and width of the current 


system as 


and 


Tus, thelrotealh current d is 


These definitions for width and J are important, especially 
for model current systems containing ionospheric sheet 
currents oriented along the auroral oval, in which case 
neither the width or latitudinal current distribution are 
necessarily constant with respect to longitude. 


Writing equation 2.6 in differential form 


dd = (0) )r, dO, (93 
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(dJ is the current flowing in the primary loop located at 
6.) and using eq. 2.3 and 2.6, the total contribution to 


each field component can be expressed as 


S 
BilrgByetg) =k f 1(6,) | f ar r, de, 
Oy Toop(6,) 
ag 
where is the line integral over the primary loop 


loop(6, ) 

located at ory 
To evaluate the line integral in eq. 2.8 the ds 

element at any point on the loop must be known. Therefore, 
to accomplish this, we separate the E-W loop into elements 
a,b ,c and d (see Fig. 2.2a) and derive an expression 
for ds on each element. First we derive an expression for 
ds along the field line passing through point P(r, 56, 5,) 
in the ionosphere. This expression is then substituted into 
BG. 2.0107 dc , thus enabling us to obtain the contribution 


to the line integral in eq. 2.8 due to field-aligned 


elements a and c 
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The equation for a field line through P(r, 58, 5) 


in the ionosphere is 


r 


k ors 
eee es ie Re he hy (2.9) 
sin’ @ 


k 


One method of obtaining the ds element along this parti- 


cular field line is by setting 





ds = |ds| i, 
= ds iy (2.10) 
where iy is the unit vector parallel to the field line and 
the magnitude | ds | is given by the equation 
4s 
ds = ‘ (3 cos” 6 +1) de (25 tg 
sin 6 ‘ 
(ds can be derived using eq. 2.9). The unit vector i 
can be found from the relationship 
oF i, x $¢ te oe) 
where 
oe, 
$ : VY (sin 6/r) (2.13) 
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THUS 


‘Fes 1 (2 cos 6 r + sin 6 8) (2.14) 


(3 cos? 6 +1) 


Substituting eq. 2.11 and 2.14 into eq. 2.10 we find 


ds = ds, r + ds. 6 + ds. 4 
where 
ds, >» 2% cin @ de 
ds, = yr de (2A 
and ds, =" 0 


Since we have expressions for ds, j ds, and ds 


dc (see eq. 2.2) is completely specified and hence, the 
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where the minus sign and 96 = dr pertain to element a , 
plus sign and 6 = bw to element c,and r (see eq. 2.9) 
and elements of matrix A (see eq. Al.12) are functions of 
6 

The solutions for ds along the ionospheric (b) 
and equatorial elements (d) are much simpler since the flow 


is along a constant latitude line in both cases. The com- 


ponents of ds are 
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where the minus sign, r = ry and 6 = Oy pertain to 
aleneni—bes—the plus Siti me = R andy “6g "12 to 
Slemented . Substituting_eq. 2.16 and. 2.18 into eo. 2.8, 
we obtain a solution for the magnetic field associated with 
an E-W three-dimensional current system as shown in Fig. 2.1. 
In the following numerical results H,D and Z 
are used to represent -B. : B., and -B, respectively. The 
numerical method we used in evaluating integrals such as 
those in eq. 2.8 is discussed in appendix A3. 
Fig. 2.4A shows the latitude profile which would 
be observed on a meridian 4° east of . for an E-W current 
system 20° long and 5° wide evaluated using eq. 2.8. Fig 
2.4B, C and D show the individual contributions to the magnetic 
field from the ionospheric, western and eastern Birkeland 


sheet currents respectively. Since the contribution from the 


equatorial sheet current for this particular system is only 
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2.4 Latitude profile of the total field associated 
with an E-W three dimensional current along with 


profiles attributed to each individual sheet 


current. 


8S 


8S 






of ed. 7 


= mene F = 
Pat +] . # 

D> } / q > 

= j 
iad f ’ ; 
% S | / i alt 
| ; = 
-) oo : a I 

~<| ioe ae o> : 
wie a 
j= a4 | : 





\ 
iMsEMa TLS 











4 i 7 
. \ 1 J # 2 ie | ae | \ | ue i 
be. (ahah an ie +7 | wore ¥-3 : / hy 
; AAPLUS Tapie \ ‘ VJ / io | a a \ j , a 
; fe ; a ate | ta yy ae 
| . | 9 - S y f er 
j +e mt : oe 
\ \" of LA; A [? i 
| = / : | od J WV“ a al 
— Ane” | ls 
}——»-—_»-__- —- + ~ hd pe ———— oe 


oe oe er oF 28 ge 2 ot e at or 23 ts at 





ee Oo a nw pln a 
ae 4 / ‘ — : 

; \ ; 
ewer @ \ 





[MLEW214A (Cents) 


7 
a rs : : = a 5 a | 
— - 
iad ————— — oe — a —— 
— - - = +o 7 
=, a — ee | ‘ 
~~ \ ” tas Oe sae ps men yet = : 
hee” ie PA : 
= a a ; 

) iss & 
vast2as 1 a ymaTes 7s 
aA TazHe \ td : THRPRD T3AaK2 , ao 

a 4 Aw. 

: - Pal } oa ; iy a 

| Pas iS ‘ -_ 

\ - 
| i “~ , 


+ 
/ } 
\ 4 ; 
: \ / i ; | : 
ig % 
ee ee —— 
avg ——e oe —_——~ s 


2 06 ES or 2a ae 22 28 08 er ar 2 da a 
QU SOUT ITAL 







: | ; 
acs 

pt donee brett Istos ai to. sfitorg ‘sbughted | “ 
ae er TOL ee ie oe 

Mn One tansiith 2 = W-2 ne de 


wine 
arss 
5 7 HF 





2 a a 
VID ¢ 36 =] OJ batud ia a8 2 ovdg 
7 


7 Inst ' 


3) 


of the order of ly at these latitudes, it is not included 
separately in this figure. The contribution to the H 
profile from the Birkeland currents can be very large as 
evidenced in Fig. 2C and D (»300y at 6.) and thus substan- 
tially reduces the effect of the -H regime due to the 
ionospheric sheet current. North and south of the current 
System the +H contribution from the Birkeland currents 
dominates over the -H contribution of the ionospheric 
electrojet, resulting in +H regimes at low and high 
latitudes. 

The D profile in Fig. 2.4A represents the differ- 
ence between two large contributions from Birkeland currents 
positioned east and west of the central meridian. Since the 
latitude profile is 4° east of >. » the field from the 
eastern sheet is the dominant one. Therefore the D profile 
is negative south of 0. and positive north of 0. y THTe 
Opposite is true for a profile west of . 

Although the main contribution to the Z profile 
is due primarily to the ionospheric sheet current, a 
significant contribution is caused by Birkeland currents 
(x70y at 6.)- The effects of Birleland currents on the Z 
profile are discussed later in this chapter in conjunction 
with E-W asymmetries in the flow of Birkeland currents into 


and out of the auroral zone. 
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Fig. 2.5 shows a series of profiles at various 
longitudes east and west of . for the same current system 
as shown in Fig. 2.4. Figs. 2.5D and 2.5E illustrate the 
N-S asymmetry introduced in the Z profile east (or west) of 
the end of the ionospheric portion of the current system. 
This proves to be useful in the analysis of magnetic obser- 
vations. In Fig. 2.5F we show an example of a profile west 
of ds to illustrate the reversal of the sign of the D 
profile. 

Contour plots of the H , D and Z components as 
a function of latitude and longitude are shown in Fig. 2.6. 
We have found that such plots as these are very useful for 
mapping the positive and negative regimes of each component 
for various current systems. They are also useful in allowing 
easy visual comparison of the magnetic perturbation pattern 
of theoretical current systems with and without Earth induc- 


tion effects taken into account. 


2.2.3 North-South (N-S) current system 


A model of a N-S three-dimensional current system 
is illustrated in-Fig...2.7->> "In thiS system the Birkeland 
currents flow from the magnetospheric equatorial plane to 


the northern edge of the southward directed ionospheric sheet 
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2.5 Series of latitude profiles at various longitudes 


east and west of the central meridian. 
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Fig. 2.6A Contour plot of the H component for an E-W current 
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current and out at the southern edge to the equatorial plane 
where the path of current flow is completed by a sheet 
current directed away from the earth. We can obviously 
reverse the direction of current flow to give a northward 
directed ionospheric sheet current. 

Since the mathematical development we presented 
for calculating the field due to an E-W current system is 
Similar to that required for a N-S system, we will only 
present a brief discussion concerning the latter system. 

The ionospheric element of the N-S primary loop 
(see Fig. 2.2b) we use in the evaluation of the field due 
to the model shown in Fig. 2.7 is presented in Fig. 2.8. 

In this case we obtain an analytic expression for the field 
of this primary loop as a function of d, . Hence, the 

total field due to an entire N-S current system can be 
evaluated by integrating this expression over the length 

bw to op . The definitions of length and width of the 

N-S current system are the same as those for the E-W system. 
We also define an expression similar to eq. 3.6 for the total 
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where I(¢,) is the longitudinal current distribution 
function. Therefore we can write an expression for the 


total field similar to eq. 2.8 as follows: 


OF 
By (1528, 9%) = Kf I($,) ef d de, ry sin TF do, 
e Toop(,) 
Caray 
where f is the line integral over the N-S 
loop(¢,) 


primary loop located at op 

Since we have already calculated the line integral 
over an E-W primary loop (see eq. 2.16 and 2.18), we do not 
believe that a similar calculation for the N-S primary loop 
is warranted here. 

Fig. 2.9 shows the latitude profile 2° of %% 
for a N-S current system re long and 4° wide. Again we show 
the individual contributions to the magnetic field from the 
separate portions of the current system. Fig. 2.10 shows 
latitude profiles east and west of do i) LiesRenroriie 15 


symmetric in longitude with respect to do. while the H and 


Z profiles change sign from one side of the central meridian 
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to the other. The D component for the N-S current system 
is most interesting since the magnetic perturbations asso- 
ciated with eastward and westward travelling loops and surges 
are of a similar nature. This similarity is discussed further 
in Chapters V and VI. 

Contour plots of the three separate components 
are shown in Fig. 2.11. Although this particular N-S 
current system is probably much less in longitudinal extent 
(length) than those observed in the ionosphere and magneto- 
sphere, these particular contour plots show the important 
aspects of any such current system with finite length. For 
example, from Fig. 2.10C and D we see that the Z and H com- 
ponents are the dominant ones to the east and west of the 
ends of the ionospheric sheet current. For a much longer 
N-S current system with a length of 60° (v2550 km), width of 
2 (v450 km) and a total integrated current of 10° amps 
(corresponding to a sheet current density of 0.4 amp/m), 
the maximum value of Z at 30° east of d. is -106y com- 
pared with .13y for the maximum value of D. 

Another interesting feature is the depression of 
the D component east and west of the current system. These 
negative D regions are due to the fact that the D component 


field from the Birkeland sheet currents dominate over the 
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~N-S CURRENT SYSTEM — 
H-COMPONENT 
N-S EXTENT 65.5°—69.5° 


E-w EXTENT 4° 
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Fig. 2.11A Contour plot of the H component for a N-S 
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this contour and the one shown in Fig. 2.6B. 
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Fig. 2.11B Contour plot of the D component for a N-S 
current system. 
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positive D field from the ionospheric sheet current (see Fig. 


2, Brey andy gD). 


2.3 Induced currents in a conducting Earth 


Since the variations of the Earth's magnetic field 
contain contributions attributed to induced currents in a 
conducting Earth as well as those from external currents, any 
model current system necessary to explain these observations 
Should contain an induced current system in addition to an 
external current system. Although the conductivity structure 
of the Earth is extremely complex (Rikitake, 1966), we must 
assume a very simple model for this conductivity structure 
in order to make the mathematics tractable. 

Bostrom (1971) used the image dipole method to 
obtain the field due to currents induced by the three- 
dimensional current system discussed in Appendix Al.2. He 
represented the conductivity structure of the Earth by a 
spherical shell of an infinitely conducting medium of radius 
bond bak sone the radius» of, the, Eanth). Thesresultsmf this 
solution are presented in Appendix A2.4. 

Ashour (1971) derived a set of relationships 
between the field produced by currents induced in a uniformly 


finite conducting spherical shell and the field of an arbi- 
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trary external current system provided the external field 

is known numerically. We have presented a much simpler 
derivation of Ashour's results in Appendix A2.2 for the case 
of infinite conductivity. The relationships between the 
internal (induced) and external (inducing) components of the 


magnetic field (equations A2.18, A2.30 and A2.31 in appendix 


A2.2) are 
B (as0-30, 7 = lypbplege (r_ 6 ) (eurag 
] ao a a ] Ble 0°%o : 
a 
0 
AUeMER Moov ety the eybernga hor tropes) comppnenfs nor; 
EE ea.) a ee Seu 6 ab ) 
) 
e 
. Bo 3 (r 28529, )dr, O2:pB2) 


where (2,65) represents the observation position on the 


surface of the Earth (r.= a) and 


b being the radius of the spherical shell (see Fig. A2.1) 
in Appendix A2.2). 
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The vertical component of the induced field given 
by eq. 2.21 can be evaluated by evaluating the vertical com- 


ponent of the external field at 2 Thiseas at soctrves Tor 
the horizontal components with regard to the first term in 

eq. 2.22. We see however, that the evaluation of the second 
term in eq. 2.22 is more difficult because of the integration 
with respect to Char If an analytical expression for the 
external field exists such that the integration can be per- 
formed analytically, then eq. 2.22 becomes readily applicable. 
If this integration cannot be accomplished analytically, 
numerical methods would then have to be used. Since the 
evaluation of the external horizontal components normally 
requires numerical integration, the additional integration 
over or. would probably necessitate the use of large amounts 
of computer time. We have found that the horizontal compo- 
nents of the external field given by eq. 2.1 can be inte- 
grated analytically with respect to Atta thus allowing us 

to obtain a solution for the induced field due to an arbitrary 
Current system external to the conducting shell. The complete 
solution for the three components of the induced field can 


be expressed as 


: 3 : - 
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(e,=-1, Eo 371) C223) 
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i] it ~ | 


where aoe and des are elements of matrices dc and 


dc respectively and are defined in Appendix A2.3 (see 


ba dim 


Appendix A2.3 for the derivation of eq. 2.23). Since most 
of the variables in dc are also common to dc A. Tt 
the inclusion of induction effects in the model current system 
discussed in the previous section increases the computer time 
required by less than a factor of 1/10 

The induced field and total field (sum of both the 
induced and external fields) for various depths (a-b) of 
the superconductive shell below the surface of the Earth are 
presented in Fig. 2.12. The physical parameters of the E-W 
external current system are the same as those for the system 
described in section 2.2.2. The latitude profiles in Fig. 
2.12 pertain to a meridian 5° east of 6 and thus can be 
compared with the profile (representing the field due to the 
external system) shown in Fig. 2.5B. The main feature 
readily apparent in this series of profiles is the large 
value of the induced Z component at both high and low latitudes. 
As a result of this effect the total Z component is highly 
attenuated at these latitudes. This effect was reported 


previously by Bostr6ém (1971). 
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Fig. 2.12 Total field (sum of both the induced and external 
fields) and induced field for various depths of 
the superconductive shell below the surface of 
the Earth computed for an E-W current loop. 
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Although the horizontal components of the internal 
and external fields are normally of the same signyethis is 
not, however, true in the transition regions from negative 
to positive values of these components. This becomes more 
evident the greater the depth of the conductive layer. 
Listed in Table 2.1 are the ratios of the magnitude of induced 
field to external field at various latitudes for the three 
Cases suown iff Fig. 2+le. Ihe interesting feature apparent 
in this table is the large variation in these ratios with 
respect to latitude, the smaller ratios occuring near the 
latitude of the current system. 

To illustrate the effects of induction as a function 
of longitude we present a series of latitude profiles at 
various longitudes east and west of d. tm Fig. 02327 and 
contour plots of the three components in Fig. 2.14. These 
figures correspond to Fig. 2.5 and 2.6 which show similar 
contour plots where induction is not taken into account. It 
is evident from these figures that the Z component is highly 
attenuated away from the ends of the ionospheric portion of 
the current system as well as at high and low latitudes. 
Another interesting feature is the small effect that the 
induced field has upon the spatial extents of the positive 


and negative regimes of all three components. We also see 
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Tabel.@.—2+-1 RatoSot—carkey nal “tor externd ym © bd 


CF ep Gmop Fore (CH? Ae + 7°54 \y, 


Depth of Superconductive Layer (a-b) 


Latitude 100 km 200 km 300 km 
55° 0.836 0.691 0.560 
60 0.736 0.526 0.375 
65 0.416 0.221 0.133 
67.5 0.452 0.238 0.142 
70 0.400 0.210 Of..127 
75 0.751 0.532 0.377 
80 0.874 0.725 0.590 


84.5 0.924 0.816 0...7.09 
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Latitude profiles of the total field associated 
with an E-W three-dimensional current system (see 
Fig. 2.4 for the external field). 
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Fig. 2.14A Contour plot of the H component including the 
induced field (see Fig. 2.6A). 
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Contour plot of the D component ‘including the 
induced field (see Fig. 2.6B). 
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that the locations of the maxima and minima for each compo- 

nent remain approximately the same. This fact is important 

in our discussion in Chapter III in conjunction with the use 
of latitude profiles to infer the physical parameters of 


the current system. 


2.4 A model current system coincident with the auroral oval 


Since the westward electrojet is observed to flow 
along the auroral oval (Feldstein, 1963; Akasofu et al., 
1965), we have modified the E-W current system to include 
such current flow. Kamide and Fukushima (1970) considered 
current flow along the model auroral oval given by the 


expression 


5° = oe MP (0, : 61) I t'cos"> (2.24) 


wherein 86 and 6 represent the geomagnetic colatitudes 
at midnight (% = 0) and noon (9 = 180°) respectively. Their 
model current system was confined entirely to the ionosphere 
and thus did not contain Birkeland currents. Although we 
developed a three-dimensional current system using their 


model oval, we prefer a parabola representation of the oval 


instead. The advantages of such a representation will 
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become apparent in the ensuing discussion. 
Since we are mainly interested in modelling the 
nightside portion of the oval, the following expression is 


used: 


ae 2 (2.25) 
Cpe cos? > + ae sin o] 


where a is the midnight (¢ = 0) colatitude and b_ the 
dawn or dusk (@ = 90° or 270°) colatitude of the auroral 
Oval. Four different parabolas are illustrated in Fig. 2.15. 
The different arrows show possible paths of current flow. 

For example we may wish to have the current flowing along 
parabola #3 in both the morning and evening sectors or 
possibly along #2 in the morning sector and #4 in the evening 
sector. By utilizing configurations such as these we have 
thus greatly enhanced our ability to model observed magnetic 
perturbations. Fig. 2.16 shows a parabola (a = 22.5° and 

b = 17.5°) superimposed upon an auroral oval (Q = 2). The 
dayside portion of the oval can also be modelled using a 
parabola that maintains continuity with the nightside portion 
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Diagram of four possible parabolas to represent 
the path of current flow in conjunction with 

the auroral oval. The two sets of arrows repre- 
sent possible paths of current for a single 
current system. 
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Fig. 2.16 Diagram showing the relationship between the 
auroral oval (Q = 2; see Whalen, 1970 and 
Feldstein, 1967) and the parabola corresponding 
to a = 22.5° and b = 17.5° (see eq. 2.25). The 
shaded area represents the ionospheric sheet 
current associated with an E-W three-dimensional 
current system modified to flow along this 
particular parabola. 
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The shaded portion centered about 6 = ‘ce repre- 
sents the ionospheric portion of a modified E-W current 
System coincident with the model oval. Latitude profiles 
at various longitudes denoted by A through F are shown in 
Fig. 2.17. The important effect evident in this series 
of profiles is the deformation of the D component profile 
for profiles east and west of the central meridian. This 
deformation (southward movement of the D component crossover) 
is caused by the ionospheric sheet current flowing at some 
angle (other than 90°) with respect to the profile meridian. 
This effect has been observed with our line of magnetometers 
and in most cases can be explained in terms of currents 


flowing along the auroral oval (see Chapters IV, V, and VI). 


2.5 N-S and E-W (latitudinal and longitudinal) current 


distributions associated with the E-W model current 


system 


2.5.1 Latitudinal current distributions 


To obtain the numerical results presented in 
section 2.2.2 we represented the latitudinal current dis- 
tribution function by a constant. From preliminary analysis 


of magnetic perturbations associated with polar electrojets 
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we found that more complex distributions of current are 
necessary to explain the observed latitude profiles. There- 
fore, we modified our model current system programs to 
include various current distributions. 

The functions we found to be most useful for 


representing the latitudinal current distribution are 


ix) Mh Fy tidy aa°> (2.26) 


and 


yt Ta pend * gy isk) (Poa 


where variable x represents a change in variable 6 such 


that the limits of integration change from 8 and 6 


-] and +1 respectively for integration over the width of the 


to 


current system (see Appendix A3, eq. A3.3). 

Some examples of various distributions of currents 
utilizing @q. 2.26 and 2.27 are illustrated in Fig. 2.18. 
The first profile (Fig. 2.18A) is an example using eq. 2.26. 


The parameters for this model system are I, a nee I, = ] 
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Latitudinal profiles associated with current 
systems having various latitudinal current 
distributions. 
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? , and the profile 5° F 


a.8,.% ,_bengthds 20° » width = 5 
of the central meridian (the current distribution is always 
normalized to give the required total integrated current). 
Eor.the second profile (Eig, 2,.18B) we used eq, 2.27 with 
the following parameters: I, S20) O08 ay I, Sa) ek Se ns S, 
b = 7.0 , length = 40° , width = 10° and the profile 10° £ 


0 
of the central meridian. This particular distribution 
function can thus be used to represent either a northern 

Or southern border enhancement. Fig. 2.18C shows the 
magnetic effect associated with an eastward and westward 
electrojet, each having different central meridians. This 
System of currents is discussed in more detail in Chapter IV. 
Illustrated in Fig. 2.18D is a double westward electrojet 
obtained by adding two current systems, each having the same 
parameters as those for the system in Fig. 2.18A. Note the 
partial cancellation of the Z component in the region between 
the two electrojets. The opposite effect is readily apparent 
in Fig. 2.18C in which case the Z components are of the same 
sense and thus add to give an enhanced -Z regime between the 
two electrojets. Both compound current systems are very 


common in the auroral zone. 
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2.5.2 Longitudinal current distributions 


In section 3.2.1 we represented Birkeland currents 
associated with the E-W three-dimensional current system as 
Sheet currents located at the eastern and western ends of 
the ionospheric electrojet. Since, if such currents exist, 
they are probably distributed all along the electrojet, we 
have modelled various distributions of field-aligned current 
flow using both the method developed by Bonnevier et al. (1970) 
(see Appendix Al.3) and the Biot-Savart method presented in 
this chapter. Eq. 2.26 and 2.27 were also used to represent 
the longitudinal current distribution function. 

Fig. 2.19 shows the results obtained using the 
Bonnevier et al. method. Eq. 2.27 was used with I, ="0 7003)", 
I, = 1.0 , a= 2.0 and b= 5.0 . No induction effects were 
included in the numerical evaluations. The shape of the 
distribution function is shown in Fig. 2.19A with the location 
of the various profiles designated by arrows labelled A 
through F . The latitudinal current distribution is uniform 
but could also have been of the forms discussed previously. 
From Fig. 2.19C and D we see that the sense of D changes 
somewhere between 5° W and 0° ($6 = 0° is normally defined 
as the central meridian). The longitude at which D cancels 


completely can be determined from the moment of the current 
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distribution and in this case is » 4° W of ¢ . Hence, in 
this particular model 2.4 x 10° amp flow into the auroral 
zone along much of the eastern portion of the electrojet 

and out along the western portion (» 6° in length). We can 
also see what at first appears to be an enhanced Z component 
south of the electrojet in the longitudinal region of 
enhanced Birkeland currents. This apparent enhancement of 
the Z component is actually caused by a reduction in the -H 
regime beneath the ionospheric portion of the current system 
due to the enhanced +H regime attributed to the Birkeland 
currents (see Fig. 2.4). The possibility of having such an 
effect occur in the auroral region complicates the interpre- 
tation of magnetic perturbations associated with polar elec- 
trojets. Also evident in Fig. 2.19 is the E-W asymmetry in 
the D component (larger values of D occurring at the western 
end of the electrojet). 

In Fig. 2.20 we demonstrate how longitudinal 
distributions of Birkeland currents can be included in a 
model current system by using Biot-Savart's law. Various 
distributions of E-W current systems can be added together 
as illustrated in the diagram at the top of the figure. In 
this example induction effects are included, thus making the 
negative peak in the Z component south of the current system 


less than the maximum negative value of the H component. 
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Since the ionospheric portion of three-dimensional 
current systems is dependent on the various conductivities 
in the ionosphere, we studied the effects of variations in 
height on model current systems to ascertain whether an 
integration over height is warranted. Fig. 2.21 shows a 
series of profiles for the same current system at various 
heights above the Earth's surface. It is evident from 
these profiles that the overall shape of the profiles remains 
relatively unchanged over this height interval with only a 
variation in magnitude on the order of 15 to 20% . Thus we 
feel that such an integration is not warranted. 

It is interesting to note that the distributions 
discussed in this section are also applicable to the N-S 
three-dimensional current system. More complicated distri- 
butions such as shears in E-W and N-S current systems are 


discussed in later chapters. 
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CHAPTER III INSTRUMENTATION, DATA PROCESSING AND ANALYSIS 


3.1 Instrumentation 


During fall 1969 and summer 1970 we installed a 
line of three-component fluxgate magnetometers in western 
Canada along the corrected geomagnetic meridian of 301°E 
(see Fig. 3.1 andiTable 3.1 for the coordinates of the 
observatories and the definition of corrected geomagnetic 
coordinates). The nine observatories were closely spaced 
in the latitude range 58.7° to 77° N, and data to extend 
our line were obtained from standard observatories at 
Newport (55.3° N) and Resolute Bay (g4.3° N). The coordi- 
nates and code names of the complete line of observatories 
are given in Table 3.1. 

Since, at most, only five stations were operative 
simultaneously during fall 1969, all the data used in the 
present analysis were obtained during summer 1970. Thus the 
description of the instrumentation pertains to that used 
during the 1970 field season. 

Eight of the nine stations recorded three compo- 
nents (H, D, Z) of the magnetic field variations in digital 
form on magnetic tape while a paper chart recorder was used 
at Crow Lake (sample rate of one data point for each compo- 
nent every 9 seconds). The sample rate for the digital 


system was one data point for each component every two 
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Table 3.1 Locations and code names of the magnetometer 
Stations. 


Code Corrected Geomagnetic Coordinates* 
Observatory Name Latitude(°N) Longitude(°E) 


Resolute Bay TS 84.3 306.0 
Cambridge Bay CAMB a Oe 300.3 
Baker Lake BAKE foot 320.4 
Fort Reliance FTRL 71.4 300.0 
Fort Smith FTSM 68.1 299.8 
Fort Chipewyan maa 57 2 BA iss 
Fort McMurray F TMU aE 302.7 
Crow Lake CRLK 53.65 Ua 2 
Meanook MEAN WAR? sO] 2 
Leduc LEDU Doge et fe ae 
Calgary CALG 5B. / sUr et) 
Newport NWPT ae 299.5 


_—_:. rr 


*The corrected geomagnetic coordinate system is a non- 
orthogonal system derived using the first five terms of 
the spherical harmonic expansion of the Earth's magnetic 
field (see Hakura, 1965). 
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seconds, and timing was accurate to within + 2 sec. The 
dynamic range of the system was + 1000y with a sensitivity 
of + ly. The response was flat in the frequency range 

0 to 0.3 Hz, with the high frequency end of the spectrum 
being set by the characteristics of the aliasing filters. 

A more complete description of the instrumentation is given 
in Appendix A4. 

To complement our line of stations we installed 
all-sky cameras (on loan from the University of Alaska) at 
Meanook and Fort Smith. Photographs were taken at 1 min 
intervals with an exposure time of 8 sec on 16 mm Kodak 


Tri-X Reversal film. 


3.2 Data processing 


The process of transferring data from field tapes 
(7 track) to nine track tapes for use on the computer proved 
at times to be a long and tedious task. When the data was 
rewritten on nine track tapes the WWVB time blocks (see 
Appendix A4) were interpreted using a special WWVB read sub- 
routine. This allowed the beginning time along with the 
station code name and reference number to be included in a 
label before each block of data. Reference numbers '1970' 


and '1111' were used to designate whether the block of data 
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was complete ('1970') or incomplete ('1111'). Short or 
incomplete blocks were usually attributed to power or 
equipment failures. When two or more bad blocks occurred 
between successive WWVB blocks, the good data blocks in 
between these bad blocks could not be timed automatically 
on the computer. If it became desirable to use such data, 
an attempt was made to correlate micropulsations occurring 
in this time interval with those recorded at adjacent 
stations. In most cases this method proved to be satisfac- 
tory. When the WWVB signal was too noisy to be interpreted 
on the computer, a Cal Comp plot was made to allow for a 
possible interpretation manually. 

After the data from all stations were transferred 
to nine track station tapes, a series of event tapes was 
created with each tape containing a series of files. In 
each file all available data from the line of stations for 
a specified period of time (event) were written using the 
same format as was used for the nine track station tapes. 
This allowed the study of a specific event such as a polar 
magnetic substorm or micropulsation train by only utilizing 
one file rather than eight station tapes. 

Most of the D component magnetic data recorded at 
Meanook were in error due to a faulty D component amplifier. 


The gain of the amplifier changed erratically during the 
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summer and thus the output had to be compared with standard 
magnetograms from Meanook. Fortunately during most of the 
events studied the gain remained at a constant (larger than 


normal) value and therefore was easily corrected. 


3.3 Data analysis 


In the study of spatial and temporal development 
of polar electrojets we used five different methods of pre- 
sentation of data from our line of stations and supporting 
stations in the middle and high latitudes of the northern 
hemisphere. To complement the discussion of these methods 
examples of each can be found in Chapters IV, V and VI. In 
the ensuing chapters we used H, D and Z to represent the 


perturbation field. 


3.3.1 Component magnetogram format 


In this method the H, D or Z components of the 
magnetic data from our line of stations for a specified 
period of time are plotted with the northernmost station 
at the top of each plot. The magnetograms from supporting 
stations (NEWP, RSLT and CRLK) were hand digitized and, if 
desired, could be included in these component magnetograms. 


This method of data presentation is used mainly to give a 
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qualitative interpretation of the specific event being 


studied. 


Seon 2 Component contour diagrams 


A more quantitative method of presenting the data 
from our line of stations is accomplished by plotting equal 
intensity contours of each magnetic perturbation component 
with universal time and latitude as coordinates (see Zaitzev 
and Bostrom, 1971). The Z contour diagram is especially 
important since the approximate latitude and width of the 
electrojet can be obtained from the Z = 0 contour and the 


distance between lines of maximum and minimum Z respectively. 


3.3.3 Three-dimensional (3-D) diagrams 


Another method of presenting the same data contained 
in contour diagrams is drawing perspective views of three- 
dimensional surfaces. This is accomplished with the use of 
a perspective plotting subroutine written by H. Jesperson 
(Iowa State University) and modified by R. J. Beebe (Stanford 
University). Although there is a loss of quantitative infor- 
mation as compared with that available from a contour diagram, 
we feel that such 3-D diagrams are useful especially for com- 


plex spatial and temporal developments of current systems. A 
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good analogy would be the comparison between a map of equal 
elevation contours of a mountainous region and a photograph 


of the same region. 


3.3.4 Polar plot 


The most widely used method of displaying magnetic 
data from stations distributed both in latitude and longitude 
is that of plotting a vector representing the total horizontal 
component of the perturbation field for each station on a 
polar plot. Since there are numerous coordinate Systems which 
can be represented with a polar plot, a choice had to be made 
as to what particular coordinate system would be best suited 
for ordering magnetic data at high latitudes. Whalen (1970) 
(also see Comfort, 1970) has discussed the advantages of 
using the corrected geomagnetic (C.G.) coordinate system (see 
Hakura, 1965) over other coordinate systems (i.e. centered 
dipole, eccentric dipole and geographic coordinate systems). 
Although some difficulties could possibly arise with the 
coordinate system (see section 3.4), we have displayed our 
data in the C.G. system of coordinates for all polar plots 
in Chapters IV and V. 

Whalen (1970) has also presented a very useful 


nomograph for plotting auroral ovals in the C.G. coordinate 
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System and for converting from geographic to C.G. coordinates. 
Another valuable aspect of this nomograph is the ease with 
which the C.G. local midnight meridian (defined as the meri- 
dian which passes through the anti-subsolar point) can be 
obtained for any season of the year. Thus, rather than 
representing the azimuth of the polar plot as C.G. longitude, 
C.G. time can be used instead (C.G. time being computed on 
the basis of one hour equals 15° SF7 E.G. longitude). To allow 
for a constant format for all polar plots with the position 
of stations remaining the same, we chose to represent the 
azimuth as C.G. longitude with C.G. noon being indicated by 
a large dot which rotates around the polar plot as a function 
of universal time. The size of the polar plots are such that 
they can be used in conjunction with Whalen's nomograph. 

Fig. 3.2 shows the distribution of magnetic stations 
used in the present analysis (see Table 3.2). The C.G. 
coordinates for these stations were obtained using tables 
calculated by Hakura (1965). An example of a polar plot 
displaying actual magnetic data is shown in cTo> 2yeee rine 
position of the auroral oval (Q = 1) was determined using 
the auroral oval plotter. Although the distribution of 
available stations is rather sparse in many areas, we are 


still able to obtain valuable information in order to infer 
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CORRECTED GEOMAGNETIC COORDINATES 


Fig. 3.2 Distribution of magnetometer stations in the 


northern hemisphere that are used in the present 
study. The outline of the continents is in 


corrected geomagnetic coordinates after Whalen (1970). 
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Lap 6.3.2 Locations and code names of magnetometer 


stations used in polar diagrams. 


Observatory 


Alert 

Resolute Bay 
Mould Bay 
Cambridge Bay 
Baker Lake 
Fort Reliance 
Fort Churchill 
Barrow 
Narsarssuak 
Great Whale River 
Fort Smith 
Lynn Lake 

Fort Chipewyan 
Leirvogur 
Tixie 

Fort McMurray 
College 

Kiruna 

Crow Lake 
Meanook 

Leduc 

Sitka 

Ottawa 

Lerwick 
Calgary 

St. John's 
Lovo 

Newport 
Victoria 
Fredericksburg 
Boulder 

Dallas 

Tucson 

San Juan* 
Honolulu* 
Guam* 
Huancayo* 


Code 
Name 


ALER 
RSLT 
MOUL 
CAMB 
BAKE 
FTRL 
CHUR 
BARR 
NARS 
WHAL 
FTSM 
LYNN 
FTCH 
LEIR 
TIXI 
FTMU 
COLL 
KIRU 
CRLK 
MEAN 
LEDU 
SITK 
OTTA 
LERW 
CALG 
JOHN 
LOVO 
NEWP 
VICT 
FRED 
BOUL 
DALL 
TUCS 
JAUN 
HOLU 
GUAM 
HUAN 


Corrected Geomagnetic Coordinates 


Latitude (°N) Longitude(°E) 
86.5 122.6 
84.3 306.0 
80.5 264.3 
ad KY S00. 
ne | 320.4 
71.4 300.0 
70.0 a2, 0 
69.7 24740 
69.0 45.6 
68.2 54.0 
68.2 299.8 
Ove 3 346.1 
67.0 a04 3 
66.3 eau 
65.6 194.9 
65 802.7 
64.9 260.3 
64.3 104.8 
63.8 3014 
62.5 30142 
Glee 301.5 
59.8 2466 
58.9 $55.7 
58.9 84.3 
5847 302.0 
Saad 29.8 
55.9 97.9 
bE AG 299.5 
53.4.9 292.6 
51.8 So2ve 
49.3 21 56¢x 
43.4 Seocs 
39274 JAA 
29.6 So 
Pilea 266.5 

4.0 ote. o 
-0.6 ae ie 


*Geomagnetic Coordinates 
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some of the necessary parameters for determining an equivalent 
three-dimensional current system necessary to explain the 
observed perturbations. For example when we compare these 
disturbance vectors with those shown in Fig. 3.4 for a three- 


6 amp), we are thus able 


dimensional E-W current system (10 
to determine the approximate latitudinal and longitudinal 
position of the current system. This is discussed in more 


detail in Chapter VI. 


3.0.% Latitude profile 


From the results presented in Chapter II on model 
Current systems, it is readily apparent that the use of 
latitude profiles provides an excellent technique for dis- 
playing magnetic fields associated with those systems. Since 
our line of stations is positioned approximately along the 
a0 Teun, Ge meridian, displaying the magnetic data in the form 
of latitude profiles allowed us to make comparisons between 
observed profiles and theoretical profiles. This technique 
was used previously by Walker (1964), Bonnevier et al. (1970) 
and Kisabeth and Rostoker (1971) and has been discussed in 
detail by Kisabeth and Rostoker in conjunction with inferring 
three-dimensional current systems from data recorded with our 


line of stations. 
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Throughout the discussion of the spatial and tem- 
poral development of current systems, some characteristics 
of latitude profiles are used to infer the physical parameters 
of these systems. In the following discussion the physical 
parameters of interest (with reference to E-W three-dimensional 
Current systems) are the moment, width and latitudinal distri- 
bution of current. In Chapter II (section 2.3) we showed that 
the effects of induction do not seriously alter these para- 
meters (also see Bostrém, 1969). All of the figures referred 
to in this discussion are found in Chapter II. 

The moment of the current system is defined as the 
latitude at which the Z component is zero for meridional 
Station lines that lie under the ionospheric portion of the 
current system (see Fig. 2.5A and 2.18B). The moment cannot 
always be easily defined, however, particularly when the Z 
profile becomes distorted due to two or more maxima in the 
latitudinal current distribution (see Fig. 2.18D). ARotheY 
physical situation that will affect the Z = 0 crossover is 
the longitudinal distance between the line of stations and 
the end of the ionospheric segment of the three-dimensional 
current system. If the line of stations should lie to the 
east or west of the ionospheric portion of the current system, 


the Z= 0 crossover will be located slightly south of the 
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actual moment of the current system (see iG de «alla dasa 
Fig. 2.6C). The location of the peak in H (negative for 
westward and positive for eastward electrojets) should 
coincide with the moment of the current system. Since this 
peak can be adjusted easily and still fit the observed data, 
it does not serve as a very accurate indicator of the location 
of the moment of the current system, whereas the Z = 0 
crossover is usually well defined by our observed data. 

The width of the-current system is approximated by 
the separation between the maximum and minimum in the Z 
profile. Again this is true only for single current systems. 
The width measurements should be representative of the actual 
Width of the current system if the line of stations is near 
or under the ionospheric segment of the current system. If 
the profile is obtained for a current system far to the east 
or west of the line of stations, the separation of the peaks 
in the Z profile will be greater than the width of the actual 
current system (see Fig. 2.6C). One important factor to look 
for in the Z profile is the shape of the Z peaks. If the line 
of stations is located far from either end of the ionospheric 
segment of the current system, both peaks will be very broad 
and thus the maximum will not be well defined (see Fig. 2.5E). 


The shape of these peaks can also be affected by the latitu- 
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Cia otethroucion O1 Current aS 1) lUuStrated 1f FIG. 2.1oB. 
Qualitatively, then, the shape of the peaks may give some 
indication of the distance of the line of stations from the 
end of the current system. Also the shape of the peaks (along 
with the H profile) may indicate asymmetries in the latitu- 
dinal current distribution. 

The width of the negative H regime can also be used 
to infer the relative width of the current system; however, 
the inferred width is also dependent on the distance of the 
line of stations from the end of the current system (see Fig. 
2.6A). Therefore we must be cautious in using the width of 
the negative H regime as a quantitative measure of the width 
of the ionospheric electrojet. 

The D perturbation is positive to the south and 
negative to the north of the ionospheric portion of the three- 
dimensional current system for latitude profiles west of the 
central meridian. For profiles east of the central meridian 
the sign of the D perturbation is just the opposite of that 
discussed above (see Fig. 2.5B and F). Along the central 
meridian the D perturbation from field-aligned currents will 
be cancelled. Therefore we can use the D profile to determine 
whether the station line lies along the central meridian of 


the ionospheric portion of the current system. The D profile 
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can become highly distorted at times due to the presence of 
eastward and westward travelling surges and loops (see 
Chapter V). 

In summary the Z profile appears to be the most 
important of the three profiles for inferring the physical 
parameters of the current system or systems responsible for 
the observed magnetic perturbations, although all three are 
necessary to infer the most probable model. Also, all three 
components are necessary to study the effects of induction 
that are included in the observed magnetic perturbations (see 
Chapter II, section 2.3 and Chapter VI). The D component 
appears to be the most important for detecting disruptions 
of the auroral arc systems and also flow of current along 
the auroral oval (see Chapter II, section 2.4). 

The longitudinal distribution of current (and length 
of the current system) is more difficult to ascertain from 
observations of magnetic perturbations. A discussion of this 
problem is included in Chapter VI. 

The resolution of our line of stations is very 
difficult to determine since the interpolation between two 
stations for one component may be accomplished with the aid 
of all components north and south of this particular regime. 


A discussion of this problem is presented in Appendix A5. 
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Since the arc spacings are of the order of 30 to 40 km 
(Akasofu et al., 1966) with arc thicknesses varying from a 
few hundred meters to 10 km with an average thickness of 1 
km (Magge and Davis, 1968), our line of stations definitely 
does not have the necessary spacing to resolve individual 
arcs within a group of arcs having this average spatial 


distribution. 


3.4 Baselines and the orientation of HDZ coordinates in the 


C.G. coordinate system 


Since we cannot, at the present time, define an 
absolute baseline from which to measure magnetic perturba- 
tions associated with current systems, we set the baseline 
(sometimes referred to as the DC level) at each station just 
before the onset of the particular disturbance we are 
interested in analysing. To determine the DC level we 
usually used magnetograms from observatories positioned both 
east and west of our line of stations to allow for possible 
spatial effects such as westward travelling surges. We also 
used this method of setting the DC level to examine pertur- 
bations of existing current systems. Latitude profiles with 
such baselines are referred to in the text as differential 


profiles and prove to be invaluable in the study of magnetic 
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disturbances due to the passage of auroral forms such as 
travelling loops and surges. 

Since the magnetometers along our line of stations 
are oriented in HDZ coordinates (H being the horizontal com- 
ponent of the field in the local magnetic meridion plane), we 
related this frame of reference at each station with lines of 
constant C.G. latitude and longitude. This is necessary since 
we wanted to compare the orientation of current systems with 
that of auroral ovals plotted in C.G. coordinates (Whalen, 
1970). Comfort (1970) represented lines of constant C.G. 
latitude by a family of concentric ellipses (eccentricity = 
0.5914) centered about a point on the Earth's surface with 
geographic coordinates 80°N and 280°F. To relate the orien- 
tation of the magnetometer at each station to the C.G. 
coordinate system we defined a C.G. declination at each 
station as the angle between the direction of the North 
geographic pole and the normal to the constant C.G. latitude 
line. The C.G. declination was determined by calculating 
the gradient (in spherical coordinates) of the function given by 
Comfort (1970) to represent the ellipses. Table 3.3 shows 
how closely the values for magnetic declination and C.G. 
declination correspond along our line of stations. Thus we 


should be able to compare the orientation of the current system 
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Table 3.3 Comparison of the magnetic declination and 
corrected magnetic declination (defined in 
text) for olir Tine of stations. 


Corrected 
Magnetic Magnetic 
Declination Declination Difference 
Observatory (°F) (PE) Py 
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as measured in the HDZ system with the auroral oval in C.G. 
coordinates. 

Since constant C.G. latitude lines can be represented 
by ellipses on the surface of the Earth, we can expect to have 
variations in lengths of arc measured in degrees of C.G. 
latitude. This variation is dependent upon the location at 
which the measurement is made due to the orientation of the 
major and minor axes of this family of ellipses. Such errors 
could be serious when using the latitude profile techniques 
to infer actual widths of current systems. We were fortunate 
to have our line of stations located in a C.G. longitudinal 
regime such that this error is negligible. This would not, 
however, have been the case had our line of stations been 


positioned in Eastern Canada or in the Alaskan sector. 
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CHAPTER IV DYNAMIC DEVELOPMENT OF THE POLAR ELECTROJETS 


Ae  inerodauction 


The dynamic development of two polar magnetic substorms 
are included in this chapter. Before we present these two sub- 
storms, the results of our study of the eastward electrojet in 
the evening sector are discussed with emphasis on the response 
of this electrojet to substorm activity in the midnight sector. 
Heppner (1954) has claimed that positive H component bays in the 
evening sector are associated with the development of an east- 
ward electrojet at auroral zone latitudes, while Akasofu et al. 
(1965) contended that these positive H bays are caused by the 
return current from a westward electrojet flowing far to the 
north. We hope that the results presented in this chapter will 
help to cast some additional light on this controversy. 

In section 4.2 the magnetic field due to a combina- 
tion model current system including both an eastward and west- 
ward electrojet in the auroral zone is discussed in detail. 

We feel this presentation is necessary to better understand the 
magnetic perturbations observed in the evening sector as dis- 
cussed in section 4.3. In sections 4.4 and 4.5 the dynamic 
development of the substorm events that occurred on June 15, 


1970 and July 14, 1970 respectively are presented. 
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4.2 Magnetic field associated with a model three-dimensional 


current system containing both an eastward and westward 


electrojet 


To better understand the magnetic observations 
discussed in the following sections of this chapter, we pre- 
sent a further discussion of magnetic fields associated with 
a model current system containing both an eastward and west- 
ward electrojet as shown in Fig. 4.1A (see Chapter II, section 
2.5). Represented in this figure are the ionospheric portions 
of two three dimensional E-W current systems with oppositely 
flowing current. The central meridian of this combination 
Current system is represented by the 0° meridian. The center 
(6. and ) of the eastward current system is (62.5°, i tah 
while the center of the westward current system is located 
to the north and east (67.5°, +15°) of the eastward system 
(see Chapter II, section 2.2). Both electrojets are 40° 
in length and 5° in width. The latitudinal current distri- 
butions shown in Fig. 4.1B, C and D for each system are of 
the form given by eq. 2.26 in Chapter II. The effects of 
induction weresaligsontncluded by assuming a superconductive 
layer at a depth of 250 km. The total integrated current in 


6 6 


the system is 2 x 10° amp (10° amp for each electrojet). 
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The latitude profiles shown in Fig. 4.1B, C and 
D are for meridians = Ee hg and +15° respectively. In 
Fig. 4.1B the effect of the westward electrojet appears as 
a slight enhancement of the -H regime north of the eastward 
electrojet. Since this profile is coincident with the 
central meridian of the eastward electrojet, the contribution 
to the D profile is solely from the Birkeland currents 
associated with the current system containing the westward 
electrojet and vice versa in the vicinity of the central 
meridian of the westward electrojet (see Fig. 4.1D). The 
Z component in Fig. 4.1B is almost completely due to the 
eastward electrojet, although the small +Z regime to the 
north of the -Z regime is due to the +Z regime off the end 
of the westward electrojet (see Chapter II, Fig. 2.14C). 

Fig. 4.1C shows the profile ontthe 0° meridian of 
the overall system. In this example the -H and +H regimes 
resemble one another very closely. The Z profile shows a 
pronounced enhancement between the two electrojets due to 
the addition of the -Z fields from both electrojets. The D 
profile is greatly enhanced (both the negative and positive 
regimes) and is caused by the fact that the Birkeland currents 
associated with both current systems near this particular 


meridian are flowing out of the ionosphere (see Fig. 4.1A). 
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THEGRETICRL PROF Je 
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Model three-dimensional current system containing 
an eastward and westward electrojet. The latitude 
profiles show the magnetic field along the £15 2e¢ 
0° and +15° meridian. 
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In the profile shown in Fig. 4.1D the field from 
the westward electrojet dominates both the H and Z profiles 
as was the case for the eastward electrojet along the Sy nite 
meridian. As mentioned previously, the D profile is due 
entirely to the Birkeland currents associated with the east- 
ward electrojet:', however, in the longitudinal regimes of 
-35° and +35° the Birkeland currents located in these areas 
determine the relative shape of the D profiles. 

Fig. 4.2 shows a polar plot of the magnetic 
perturbation vectors associated with the current system 
illustrated in Fig. 4.1A.° This figure clearly depicts the 
enhanced +D and -D regimes in the region where the two 
electrojets overlap in longitude (the 270° meridian in this 
plot corresponds to the 0° meridian in Fig. 4.1A). If we 
assume that 270° represents the midnight meridian, then the 
total current system produces a substantial dawn to dusk 
magnetic field perturbation on the nightside of the polar 
cap. Also note the almost complete cancellation of the 
horizontal perturbation field on the dayside (see Chapter III, 
Fig. 3.4 for the distribution of horizontal field vectors 


associated with a single E-W three-dimensional current system). 
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HORIZONTAL MAGNETIC PERTURBATION VECTORS . DEPTH OF SUPERCONDUCTIVE LAYER - 250KM 
90 
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CORRECTED GEOMAGNETIC COORDINATES 


Fig. 4.2 Polar plot showing the distribution of horizontal 
magnetic perturbation vectors due to the current 
system orn Fitg J &@ Pp! 
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4.3 Eastward electrojet and associated substorm events 


of July 5, 1970 (Day 186) 


The development of an eastward electrojet in the 
evening sector on July 5, 1970 afforded us the opportunity 
to study the response of this particular electrojet to 
polar magnetic substorm activity. The H, D and Z component 
magnetograms and the associated Pi 2 micropulsations are 
Shown in Fig. 4.3 and 4.4 respectively. The Pi 2 micropul- 
Sations presented throughout this thesis were obtained by 
numerically filtering the magnetic data from our line of 
stations using a Butterworth filter subroutine written by 
Alpaslan (1968). The component magnetograms representing 
data from eighteen observatories in the northern hemisphere 
are shown in Fig. 4.5. The letters on the various magneto- 
grams allow each trace to be followed. The + and - signs 
on the Z component magnetogram for SITK were necessary since 
this particular component was accidentally reversed during 
the plotting of this diagram. Also note the change in 
amplitude scale as indicated on the right hand side of this 
diagram. 

The purpose of the following discussion is not to 
give a detailed dynamic development of the eastward electro- 


jet; but rather to show the response of the eastward electro- 
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ig. 4.3 Magnetograms showing the development of an east- 
ward electrojet in the evening sector on July 5, 
1970. 
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DAY 186 


Z COMPONENT 


0200 0230 0300 0330 0400 0430 o0soo 0530 0600 


CAMB 
ie FTRL 
Lj FTSM 
ae FTCH 
a etm 12" 
es MEAN 
Ga LEDU 
a CALG 





0200 0230 0300 0330 eh tele) 0430 osoo 0530 0600 


Teor 
al: 


BOL OMPONE NT 


(200 0230 0300 0330 nuo0 0430 0SO00 0S30 0600 
UT 


ig, 4.4 Py 2 micropulsation activity (40-150 sec) associated 


with the event shown in Fig... @.3. 
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jet to substorm activity in the midnight sector. Although 
the magnetic recording station at FTRL was not in operation 
at the time of this particular event, BAKE (75.1°) proved 
to be invaluable for the interpretation of magnetic pertur- 
bations between FTSM and CAMB. The Kp indices for this day 
Were tony 87, 4=5 ee 38, SOR e+e 3-3; ZKp = 23) and for the 
previous 9 hours (Day 185) were lo or less. 

The gradual development of the eastward electrojet 
Shown in Fig. 4.3 began at » 0245 UT, nearly coincident with 
the development of a large -D regime observed at high lati- 
tude stations (most pronounced at ALER, RSLT, MOUL and BAKE). 
The eastward electrojet as observed with our line of stations 
was very broad at 0300 UT with no well developed peaks in the 
Z profile. The latitude at which the Z component changed 
sign was ~ 66.0°. The maximum H perturbation was »~ 60y with 
the +H regime extending to low latitudes. This is clearly 
evident in Fig. 4.3. In fact, similar +H perturbations were 
observed at HUAN, GUAM and HOLU. The development of this 
broad +H regime coupled with the large -D field observed at 
high latitudes is consistent with the pattern of the DP 2 
current system (Nishida:and Kokubunjet971). ‘Itsisrinteresting 
to note also that KIRU recorded a -H perturbation while LOVO 


recorded a very small transition H (+H to -H). This may be 
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indicative of the current cell in the dawn sector associated 
with the DP 2 system. 

Immediately after 0300 UT the H profiles from our 
line of stations showed the development of a well defined 
eastward electrojet. At 0310 UT the electrojet showed a 
maximum H perturbation of +152y at FTSM. A meaningful 
width determination of the electrojet at this time was not 
possible. By 0322 UT (see Fig. 4.6A) the eastward electrojet 
was well defined with a maximum peak of 160y and a width 
OF TO: Pa, Also, by this time the +H perturbation field at 
lower latitudes had decayed substantially. Coincident with 
this decay the -D perturbation at higher latitudes subsided 
with the large -300y" bay at MOUL having almost completely 
disappeared by ~ 0330 UT. As the low latitude +H field 
disappeared, the southern (+Z) peak in the Z profile became 
more well defined. This is what we would expect to observe 
as was shown in our model studies presented in Chapter II. 
For example in Fig. 2.18B we showed how the Z component peak 
at the equatorward side of the electrojet is not well defined 
due to the particular latitudinal current distribution used. 
If the distribution were to gradually change to the form shown 
in Fig. 2.18A, the Z peak would experience an enhancement 


along with becoming more well defined. Thus, we feel that 
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the observed profiles can be explained by the decay of the 
DP 2 system along with an enhancement of the eastward electro- 
jet in the auroral zone. This is clearly shown in Fig. 4.6A 
and B. At ~ 0419 UT the eastward electrojet was 1 6-8° wide 
with a maximum H perturbation of « -200y. 

In the sequence of latitude profiles shown in > ae 
4.6, the development of a westward electrojet is evident to 
the north of the pre-existing eastward electrojet. This 
development, in fact, was associated with the onset of sub- 
storm activity in the late evening and midnight sectors. 
The onset of a substorm just prior to 0400 UT is evident 
from the H component magnetograms for WHAL and CHUR. The 
increased micropulsation activity at this tumeris also indi- 
cative of a substorm onset. To determine the magnetic 
response to this substorm activity as observed in the evening 
sector, we set the DC level at 0357 UT and analyzed the 
ensuing development. Fig. 4.7 shows the differential profile 
and differential polar plot for »~ 0416 UT. (The length of 
the perturbation vectors in the polar plot is measured to 
the base of the arrow head.) To illustrate the approximate 
Orientation of our line of stations with respect to the 
auroral oval, a Q = 2 oval was drawn on the Prot.) Sir) 1a OB A 


and B the eastward electrojet cannot be explained in terms of 
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The development of both electrojets is depicted here. 


Note the enhancement of the D component. 


See 


Appendix A5 for a discussion of the technique of 


TITting.curves toa the data points. 







Oo 


340 


KIEFS IWLEW21LA (CUMNBE) 


a aoe 
-lETErD JALEWS1IA 1c 


0 


+ 


sdutl Ta 





109 


"4O0}ZO9S YUubLUpLW 
9yz UL WAO}Sqns e& YALM pazyeLoosse uo.das HuLuaaa 94} UL SJOfOUZDALS OMA |YY JO 


dnpLing A424 210N “IN 9L¥O “os ([ BAO Z = 0) ¥0,d ue_tod pue a, tyoud LelLQUaUasSIG Ly “Bly 


S3LUNTOYOOI ITLINOKWO39 0319394v0I 








ay ms 
L08 
SONLI Lb) _~-----e cay 
Sg og SL OL Sg 09 SS f =e Ts 
e tL = 
€ — —, ~S 
7 \\ PS A eS Ss ee, \ 
¢ Ma pice ao ~< an 
f° ¥ WG N s, 
5) / va of! Kb | N“ \ 
— 7 5 \ 
a! / oo o SS ee sures Ca s 
rg oe 
S / / oe=an ~ ES \ 
/ 
2 , 6 
= / 
m lg 
—* | OZ] 
WY | 
aS QI 
wae ¥. \ 
“Sengesne \ 
— \ 
Cy \ 
D \ 
= \ 
= \ 
WY \ 
rs) 
°93S2S ‘NIWSI YH h  9l ABO ne Wh 
Ov:9S£0 -DQd 
OsZ2 Sa y 
“33SES “NIWSIL YH h 9B8l AUN 06 0<f.2" 7 


SWIL IWSHYSAINN SYOLIZA NOT LYGYNLYSd JILINIUW TWLNOZIYOH 


-_ 
~ - -—-— < 


- ee 
- - 





(CBWWH2) 





ETEFD [MLEWUeTIA 








wv 
a OF 2a 
AOUTITAL 


T 08 





“qub{tud sft ston —.TU ATK vot (Isvo S = G) solq Slog bas stttorg Istineter 


me =. 


_ Sat nit mvotedve s dtiw betstsoe2es 10t592 pntneve sit nt etetovtzefe owe 


110 


a return current associated with a much larger westward 
electrojet flowing along the auroral Oval, especially since 
the eastward electrojet was the dominant system in the 
evening sector. In fact the existence of a westward electro- 
jet flowing along the auroral oval at 0322 UT is unlikely. 
However, in Fig. 4.7 the eastward electrojet (representing 
the perturbation of the pre-existing eastward electrojet) 
Was associated with a much larger westward electrojet (that 
had developed since » 0357 UT) flowing along the auroral 
Oval. The approximate location of the two electrojets is 
depicted by the large arrows on the polar plot (the size 

of these arrows does not represent the relative size of the 
inferred electrojets). Similar results were obtained for 
other perturbations of the eastward and westward electrojets 
in the evening sector which occurred during this event and 
events recorded on other days as well. 

The development of both electrojets in the evening 
sector for this particular event is portrayed with the use 
of the component contour diagrams presented in Pads. i cAyee 
and C. The dominance of the eastward electrojet over a 
possible westward electrojet to the north from » 0245-0430 UT 
is clearly demonstrated by the large +H regime in Fig. 5.8A. 


During this interval of time the location of the maximum +H 
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perturbation shifted slowly equatorward. This motion is also 
observable in the Z component diagram. From the slope of 

the Z = 0 contour we obtained an ‘apparent! equatorward velocity 
of ~ 150 m/sec, similar in magnitude to those obtained by 
Kamide and Fukushima (1971). The reason for emphasizing 
‘apparent’ is discussed later in this section. The rapid 
development of the westward electrojet is also shown in the 

H component diagram. From »~ 0450-0520 UT the westward 
electrojet was the dominant system as observed by our Jine 

of stations. The development of two regimes of +Z at both 
high and low latitudes and the enhancement of the -Z regime 
in between these +Z regimes after ~ 0425 UT was due to the 
development of two oppositely directed and well defined 
electrojets. The D component diagram shows the large enhance- 
ment of both the +D and -D regimes commencing just prior to 
0430 UT. This enhancement of the D component could possibly 
be attributed to :outflowing Birkeland current associated 
with these two electrojets. This idea is consistent with 

the model system discussed intthe previous section (see Fig. 
4.1C and 4.2). This effect was also reported by Kamide and 
Fukushima (1971). Associated with the large enhancement of 
the westward electrojet at ~ 0500 UT (the time of onset of 


a large -H bay at CHUR), the southern +Z perturbation de- 
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creased in amplitude and became negative even south of CALG. 
The +H perturbation associated with the eastward electrojet 
also decreased in amplitude. As the westward electrojet along 
the auroral oval decayed, both the +Z and +H regimes increased 
in amplitude and advanced poleward. From latitude profile 
analysis and the H and Z component diagrams, it is obvious 
that the strength of the eastward electrojet appeared to 
decrease significantly between 0440 and 0600 UT. However, 
the spatial and temporal changes observed may only be 
apparent changes due to the increased strength of the -H and 
-Z perturbations associated with the enhancement of the 
westward electrojet to the north. If this were the case for 
the event shown here, then this example demonstrates the 
fact that the dynamics of the eastward electrojet are not 
easily studied, especially when only a sparse distribution 
of stations is used. Also, it should be noted that the 
equatorward (or poleward) velocities obtained using the 

Z = 0 crossover motion may be meaningless when obtained 
during the development of a westward electrojet. This is 
best illustrated in Fig. 4.1B, C and D in the previous 
section. As the westward electrojet enhances (or the line 
of stations rotates beneath the current system shown in 


Fig. 4.1A, the Z = 0 crossover associated with the eastward 
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BIA to ddvoe nave ovitseen siege ud 
tetorvsoets brswtess odt Adtw bodstacaee ottadiuanahe meaty 
pnols seLavtoals brewteaw edt 2A -sbuttTqms at bezser29b oats 
beessvant zemipsy H+ bos S+ ott dtod sbeyessb Teve (sr0vus: ott 
afttovg sbutits! movi .byswelog beonsvbs bas abutelqme nt: 
ewotydo 2t t+ ,2ms%pstb tnanoqmon S bne H sit bas eteytens 

ot bevssqqs tatortosis brswiess sat to digosise edt’ tong 
eVovVeWwoH .TU 0000 bis 0820 naswisd (ltneotthngte eeseta8b 

od yfno ysm bavys2do espnsdo Tevoqmed bone -fetsaqeseas 








bns H- sit to dipnotte baessvont odd od gub gopnsdanSnerseqgs 
SAt to tnomeonsiias ois dtiw bedstooees enottedywedasqess 
40% 8269 93 s75w etdt TI . Adon ad} od tetortosls byswi2ew 
eft estsvtznomeb olaqmsxs eftdd nant . oven nwode Sneve ede 
Son 9%6 tatorvtoals byswI269 sda to 20imsnyb oat Jedd Sob 
wottudtytetb g276g2 6 vino nodw Vi fetosqze ,befbudse yl fase 
od t6id baton 9d bfuode ¢t <o2fA .bseu et enotsssecte © 
ong piteu bentstdo 2e9tsfaolsy (byswelog 10) buswrossups 


ben tstdo nodw zeetpninssm sd yem nottom. vevozeors°0°s 5 
et att ssahigs Sti biswi2ow 6 td Toemqo faved edd emtaaas ; 
2 adoivera edt mf d bas 3).af. D .pti nt besaraautit $29 
ys eedaniing tetontosts: “byswizew odt aA 





116 


electrojet moves southward even though the eastward electro- 
jet remains stationary. Therefore, we feel that such velocity 
measurements (including our previous measurement of »~ 150 m/sec) 
should be treated with care and referred to only as ‘apparent' 
velocities. If, however, the density of stations is such 

that the development of both electrojets can be ascertained, 
then (possibly through the use of modelling techniques 
described in Chapter II and VI) the true velocity and other 
important aspects of the dynamic development could be ex- 
tracted from the observed data. For an isolated electrojet 
(eastward or westward) the use of Z = 0 crossover for 

velocity calculations is probably warranted (see Cyaoter Jil, 
section 343). 

Fig. 4.9 shows some other examples of eastward 
electrojets. Fig. 4.9A and B depict the development of first 
an eastward electrojet (0504 UT profile) and then an enhance- 
ment of the westward electrojet positioned along the auroral 
oval (0603 UT profile). Fig. 4.9C and D demonstrate the 
stability of the shape and strength of a larger eastward 
electrojet associated with a large SD current system. Other 
examples of eastward electrojets are shown in Appendix A7. 

In summary, from the analysis of magnetic pertur- 


bations observed by our line of stations in the evening 
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profiles for Day 166 are in reverse order. 


Some examples of eastward electrojets. 


4.9 


rag, 


m 
is 


iy 
; 


“ie “ 


= —~833° 
elirg jeieret t 


~Sep°es 





118 


sector and by the observatories shown in Fig. 3.2, we suggest 
that the eastward electrojet has an independent component as 
well as a dependent component in relation to the polar 
magnetic substorm. Whether the mechanism for the excitation 
of both components of the eastward electrojet is the same 
requires further investigation. Due to the size of the 
dependent component of the electrojet as compared with that 
of the westward electrojet, the idea of a return current 
associated with the westward electrojet as proposed by 
Akasofu et al. (1965), cannot be ruled out at the present 


time. 


4.4 The dynamic development of the substorm event of 


June 15, 1970 (Day 166) 


In this section we present a detailed analysis of 
the substorm that occurred on June 15, 1970. The H, D and 
Z component magnetograms and the associated Pi 2 micropul- 
sations for this event are shown in Fig. 4.10 and 4.11 
respectively. The Kp indices for this day are (4-, 3+, 4, 
3+, 4=-, 3+, 3, 444 ikKple- 205).4--The -Dst-indtéees for the 
interval 0500-1200 UT are (-9, -6, -2, -1, -3, -0, 2, -2) 
indicating a lack of response of the symmetric component of 


the ring current during the time of this substorm. 
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Fig. 4.10 Magnetograms for the substorm that occurred on 
June 15, 1970 (Day 166). 
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Fig. 4.11 Pi 2 micropulsation activity (40-150 sec) associated 
with the substorm shown in Fig. 4.10. 
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From the micropulsation activity, the onset of 
the expansive phase of the substorm was identified at »~ 0702 
UT. This estimate was confirmed by the analysis of latitude 
profiles before and after the substorm onset. Before the 
onset of this particular substorm, the magnetic activity was 
relatively high as seen in the component magnetograms shown 
in Fig. 4.12 and from the high values of the Kp indices prior 
to 0700 UT. The existence of an eastward electrojet and the 
Subsequent buildup of a westward electrojet north of theteast- 
ward electrojet prior to the onset of this substorm was 
shown prevsously in Fig. 4.9. In fact a substorm was ob- 
served at CHUR around 0600 UT. Also the observation at NARS, 
WHAL and COLL of what might have been the buildup of a large 
convective type current system (SD current system), commencing 
between 0400 UT and 0500 UT, is noticeable. It is interesting 
to note that the onset of this particular substorm at 0702 
UT was hardly discernible at any of the observatories shown 
in Fig. 4.12 except for WHABL(only in the Z component). 
There does, however, appear to be a slight positive response 
of the H component at COLL associated with the onset. 

In the following analysis we feel confident that 
we were able to extract the perturbation field associated 


with the 0702 UT substorm by setting the DC level just before 
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the onset time. We are also well aware of the fact that 
the growth and decay of other systems may possibly have 
affected our results. Since the level of magnetic activity 
was so high before the onset of this substorm, evidence of 
a possible growth phase was not obtainable. 

Approximately 4 min after the onset of the substorm, 
a current system was well established with the H perturbation 
reaching -140y (see Fig. 4.13A) and the moment located at 
~ 64.0°, The current system grew steadily over the next few 
minutes while maintaining the same profile configuration, 
attaining a maximum H perturbation of ~ -200y and a width of 
~ 3.5°. Associated with this growth was a small northward 
motion of the current system (see Fig. 4.13A and B). At 
0719 UT (~ 17 min after the substorm onset) an intensifica- 
tion of the current system was observed in conjunction with 
an enhancement of the Pi 2 micropulsation activity. The 
intensification of the current system occurred on the northern 
edge of the pre-existing current system. This effeet was 
clearly apparent by 0723 UT as shown in Fig. 4.13C. Asso- 
ciated with this intensification, the moment of the resulting 
Current system moved northward » Toa By 0728 UT the ampli- 
tude of the H component perturbation under the current system 


was ~ -450y and the width of the current system was %» 4.5° 
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Latitude profiles showing the development of the westward electrojet for Day 166. 
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(see Fig. 4.13D). It is interesting to note that the D 
profile is representative of a three-dimensional E-W 
current system with a central meridian east of our line of 
stations. In Tact, up to 0729<UT the latitude Dror. 1 es. oF 
all three components are in excellent agreement with those 
presented in Chapter II. 

At 0732 UT (~ 13 min after the previous intensi- 
fication) the substorm perturbation was dramatically altered 
through the development of an independent current system at 
the northern edge of the pre-existing electrojet. The 
development of this new current system is illustrated in 
Fig. 4.14. The deformation of the profiles presented in 
this figure can best be explained by two current systems 
separated by 1 5° (Ssee-fig. 2:98D7 in Chapter \I1). The 
partial cancellation of the Z perturbation between the two 
Current systems is due to the opposite sense of the vertical 
component associated with the two current systems in this 
region. This development of a new current system caused an 
onset of what would appear at first glance to be an isolated 
substorm at CHUR (see Fig. 4.12). 

During the next few minutes the northern current 
system grew (with the H component perturbation reaching -460y) 


while the southern current system decayed substantially so 
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Same as Fig. 4.13 but showing the continued development of the electrojet. 
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that by 0734 UT the northern system was the larger of the 
two. It is interesting to note the deformation (at FTSM) 
of the.D profile exhibited in Fig. 4.14C. Such +D pertur- 
bations are indicative of E-W propagating auroral forms 
possessing N-S aligned segments such as eastward and west- 
ward travelling loops and surges. An analysis of these +D 
perturbations along with available all-sky camera photo- 
graphs is presented in the next chapter. Although no all- 
Sky camera data were available for this particular substorm, 
we believe that this development of a new current system at 
the northern border of the pre-existing electrojet was 
associated with a westward travelling surge. 

During the interval 0734-0739 UT, the two current 
systems were well defined. It should be noted that if we 
had not had a station at 67° (FTCH), the existence of two 
current systems would probably have gone: undetected. In 
fact the buildup of the northern edge of the electrojet 
commencing at ~ 0719 UT may have been due to an independent 
current system, but the resolution of our line of stations 
was not sufficient to determine whether this was the case or 
not. 

By 0704 UT, the two current systems were blended 


together forming a resultant system having a width of ~ rth 
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and a moment located at » 66.5° (see Fig. 4.15A). The 
distribution of horizontal perturbation vectors before 
0731 UT and after ~ 0744 UT showing the buildup of the second 
current system is illustrated in Fig. 4.16. The only notable 
changes in the overall pattern were the dramatic enhance- 
ments of the horizontal perturbation vectors at FTCH, FTSM 
and CHUR and a slight increase in the magnitude of these 
vectors at all other stations. The stability of the direc- 
tions of these perturbation vectors is interesting in light 
of the development of a new current at the northern border. 
There was definitely an enhancement of the eastward electro- 
jet in the afternoon and evening sectors. The horizontal 
vectors in both plots (except for the influence of the 
eastward electrojet as discussed in section 4.2) are 
remarkably similar to those shown in Fig. 3.4 for an E-W 
three-dimensional current system. The central meridian of 
the current system in these two plots is placed at 10° east 
of our line of stations. The Q = 3 oval is included in both 
plots for orientation purposes. 

The latitude profile pattern remained stable from 
0740 UT to 0748 UT after which it underwent one final altera- 
tion. Associated with a further enhancement of Pi 2 micro- 


pulsation activity, the southern border of the current system 
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drastically intensified with the maximum H perturbation 
reaching an amplitude of -870y by 0759 UT (see Fig. 4.15B 
and C). During this period of development, the width of the 
current system decreased to 1 4.5° while the moment of the 
current system moved equatorward to » 65.5° (the location 
of the moment of the current system before the initial 
intensification at the northern border). 

About 0800 UT the substorm started to decay and 
it disappeared after roughly 2 to 3 hours (see Fig. 4.15D). 
During this decay phase, the moment of the current system 
moved progressivety equatorward from 65.5° to 64.0°. 

The overall dynamic development of this particular 
Substorm as observed with our line of stations is best 
illustrated in the H and Z component contour diagrams shown 
in Fig. 4.17A and B. In the H component diagram, the pertur- 
bations of the system at 0719 UT, 0732 UT and 0748 UT are 
clearly distinguishable. The development of the new current 
system is evident from the two maxima in -H perturbation 
between 0735 UT and 0740 UT. The final intensification of 
the southern border is also well defined in this contour 
diagram. The Z component diagram shows the motions of the 
Z = 0 crossover along with the approximate width of the total 


system (designated by the dashed lines) throughout the 








fer | S 
| re ) ' oa 
nottedrws19q H mumtxem ont ddtw bart maar da 
G21.8 ett ez) TU eetO vd yO\8- to abut aS 


sat to dob tw ang .tnsigolsveb to bofr9q zat ont awd (9 bas 
edt to tHemom ait oF taw. Can vod bozso1n0b. wadeye Inara : 
nortsoof oft) as. aa. oF by6wros sups bevom satiric snort d 


[sitint oft sioted moj2ye Inavwws sdt to end to : 
.(*sbiod nyedtyon odd t6 notdsotttensdnt 
bas yso9b ot betyste myotedue sit TU 0080 two me to 


.(G21.8 .pri s92) e1wod € ot S YInpuOW yodts bovssqqs2th $f 
metey2 tnevyu9 ot to Jnomom odd .sesdq ysoobeetds patie 
°o.88 of 92.29 mor? byswrotsups vee bovom 

ict 


vefuaterag 2tdt to Fnomqolevab 3 tmanyb Tisveve saT 3 —_ 












sead ef enottss2 to seni! wo Attw bovysede 26 mrosedue 

awode emsipstb tvotn0> tnsnogmos § bns H add at betenseulit . 
-Wwiyeq sit .msypsetb tnenoqmos H 9nd al -4 brs. Khe ort at 
a16 TU SMO bas TU SEXO .TU OTTO ts mot eye edt to enotted 
tneviwo wen sit to tnemgofevab oT -sidedztientgath Da hi: 
nortsdyut1ag H~ mt smtxsm owt sit mort tnsbive etn 3 : 
Yo motseottienatnt. [snit sdT TU ODTO bas TU@EtO pie 


wotnos eda nt Sad ‘ones et > 
sheep ots Doth Jase gies s a iia 28 


4 ; 
a a 
” _ 
nan , i) 
“ f Yat | ‘ n 
if = eee ea uh Me 
‘ a we i 
, a 
—. : 7S 
7 7 — 
. oe . 
7 _ 
-_ 


a 
é. 





_ 


_ 7 
; 











ia 
a 1 


"OL6L “GL SUNP UO Pa44nd50 YeYyX WAOZSqns ay, Uoy WeUuBeLp uno UO.  UaUOdWOD H WLL'b ‘Bbq 


SWIL IWSYSAINN 
sigo o1s0 $080 0080 SSlo Oslo Sblo Ovo S£lo O€10 9210 0220 9110 O1L0 sO0lO 00L0 





Lo 


009 


seg 


os 


gl9 


ood 


99! AVG LNANOdWO9-H 


SqnLILv4 


- 





~~ \ 
iim ae 
Z, ee / 
en rr / 
 - fj j 


VY 
S\-/ 
LJ 
v ; 
A 
? 


/ 


asso ost ero =O OTO OOK 


SMIT JA@AS Vi 


<ONOT .27 snub mo beviwo20 tends miotedue oft 10? asyestb wotnos Sasnoqmoa H AND.D 








quasauddu saul, paysep aul 


fe) 


S 


fe) 


80 


0080 


——S 


sslo os2o0 


991 AVG 


SvloO 


"qal0uI9La pueMySaM Jy JO UIPLM ayy 


*yuauodwod Z ayy UOJ 3Nq YL 'p “BlL4 se awes g/lL'p 
3WIL TWSY3AINN 
Ovlo Selo o0£l0 S2Llo0 02l0 cIiZO O1loO Sol 


LNANOdWO9-Z 


aie 


te) 


5 


fo] 


OLO0 
Sl¢ 


o'o03s 


s 29 


o's 


G19 


O0L 


Gel 


osd 


o'08 


JQNLiLV7 





ite = OEtO esto «=: O80 











SS 








; MIT JAGABVIAU : 
5 _tne2e1gs 2zent! berdesb ofT .tnenogmoo § sit v0? sud AYT.& -et7 2 


+ 3 LoV7 5 


sw 


134 


development of the substorm. The major portion of the decay 
or recovery phase is not shown in these diagrams. 

In summary, the most interesting aspects of the 
development of the current system associated with this sub- 
storm are the quasi-periodic intensifications of the northern 
border of the electrojet (~ 15 min apart), including the 
dramatic development of a new current system at the northern 
border of the pre-existing system. Also the relative stability 
of the southern current system during the entire substorm is 


most noteworthy. 


4.5 The dynamic development of the substorm event of 


July 14, 1970 (Day 195) 


The polar magnetic substorm presented in this 
section occurred on July 14, 1970. The H, D and Z component 
magnetograms and Pi 2 micropulsations for this event are 
presented in Fig. 4.18 and 4.19 respectively. The component 
magnetograms for other observatories are shown in Fig. 4.20. 
The Kp indiceseTorn thisiday are. (3+, 26573", "225" 1Ts 20, 3-4 
3-; DKp = 18). The Dst indices averaged about +5 during this 
substorm with a slight decrease occurring toward the end of 


the recovery phase. 
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Tne onsetvofethis substorm was ~ 0730 UT as 
indicated by the small burst of Pi 2 micropulsation activity. 
Prior to this onset, a small eastward electrojet was observed 
centered over FIMU with an amplitude of ~ 20y (H component). 
The series of latitude profiles depicting the development of 
this substorm are shown in Fig. 4.21. By 0733 UT a well 
developed westward electrojet was observed with a maximum 
H perturbation of -90y and a moment located at 67.5° (see 
Fig. 4.21A). At this time the current system started to 
decay until (0%#45.UT at which time the northern border of this 
electrojet enhanced (1 15 min after the onset of the substorm) 
resulting in a well formed electrojet centered at » 69.0°. 
During the next 15 min the current system moved slowly 
equatorward reaching a latitude of » 68° by 0800 UT. At 
this time the electrojet was 4° wide with a maximum H 
perturbation of ~ 156y. The H component at FTMU became 
positive during this equatorward motion of the electrojet. 
The moment of the current system continued to move equator- 
ward reaching a latitude of 67.0° by 0818 UT. 

From the analysis of 10 sec resolution latitude 
profiles, the southern border of the electrojet enhanced at 
0818:50 UT (~ 34 min after the previous intensification). 


This onset is also accompanied by an enhancement of Pi 2 






7 2 ae 
- ie 
: 26 TU OEN0 » 26w myotedue “Bhd Yo teeno oAT or r- : 
-Vtivitos notfseluqotatm § 74 to fayud Fieme ont vd bedsotbat 
bevisedo 26w tetorossl9 baswi 26s TTame 6 .teeno etd. od yoryd 
. (tnanogmos H) OS » to abut if gms ie dttw UMT? are Eerenee 
to tnemgolsvsab sit pnftoraqsh eslttevq sbhustisf vorgstues ont. 
[tow 6 TU €£%0 YI .fS.8 wort mt nwone 916 mrotadue 2tay 
mumixem 6 dttw bavysedo esw tatorszels byswizow bsqoleveb 
992) °8.%3 ts bstsool tnemem (6 bis y0e- to nottsdws19q H 
ot betyste metzy2 tnsyvwo ott omit 2tdd JA .{ATS.A pra 
2fat to t9bv0d nvsdtron sit omrt dotdw ts TU 2éROC Teton Yboob 
(mvotedue sit to te2n0 ent yotts ntm af w) beonsdns tetorssefe 
"0,08 « ts bevetines totoitoste bemio? [lew 6 nf. ab ila 
vfwole bsvom moteye tnerywo sat ntm at fxn ads entiwd | 
$A =©.TU 0080 Yd Ba wn 7% obutiter 6 pnttosst bews03supe 
H mum x6 6 dtiw ebtw °b » esw ingovsueta ont omtt edd 
smso9d UMTA ts tnenoqmo> H SAT yaar wv #o nottedws1eq 
-$9(o1F3919 odd to notsom biswrotsups eft Sn iy svttt2oq. 
“10F5UPSs 9VOM ot bouNnttnoo meteye tnaviWo edt to tnemom oat 
s -TU.8180 yd °0..%2. 20 shud tier s ontdons 4 bIsBW 
shud nottufoesy og2 Of to ateylens ong mort 
2. 36 bs Bhagat t+ to a iio Me 


eee ‘ 


“ns ot a op 





i’ Ss oe 


(GAMMAS) 


FIELD INTENSITY 


(GAMMAS) 


200 


4OO 
Se ee oe 


FIELD INTENSITY 


DAY 195 


VORA) Sa MIN, 8 SEC, 





0 

















+ +--+ + + += at + + + 4 
'SS 60 65 70 75 80 8S 
LATITUDE 

DEY eS we eoMN, Se SEC, 








0 








o| 

e +—_——++ + + ——+ + + ——Hd 

SS 60 65 70 7 80 8 
LATITUDE 

Fig. 4.21 





(GAMMAS) 


FIELD INTENSITY 


(GAMMAS) 


INTENSITY 


westward electrojet. 


13y 





























a ORY 195 8 HR 21MIN. 10SEC. 
ol yi 
| 
oS 
et) 8 
T 
oO 
>t 
| 
8 | 
o——— 
ISS 60 6S 70 75 80 8S 
LATITUDE 
= DAY 195 8 HR 25MIN. O SEC. 
iB 
oO 
Lae 
<i 
Oo 
oO 
lav] 
a} 
is 
cy 
oO 
=f 
4 


oO 
aa 4+ + + 
60 





———— 


+ 
6S 70 


LATITUDE 


Latitude profiles showing the development of the 








a : : 
Ka ofaly » Fi oe CORRS ate RK 


iE ii ae 
1 me 7 Ps 2hay om 
; lov, 
ms, iw he Vs a NE 
, ee setywal 


t Nisvas: 










bagel ey 







cer ec: 
ra RTE 


= 


ELEPO TMieM2T1L 


i | 








Satay eR "a 





oe 





a 


Eis an 

CY ea og Patel ) 
30UTETAS sootirag"" id 

732 0 .wineS BH 8 zer YAO 092.2. ,KIMES AN @ = 2er 


a8 “ey 







2 ita 










S <A 
| a 
papal 
a= 
_ 


(CeWwH?} 






200 






JoMeits 





~ ®TEED 1k 


‘ meyer ih 


t-t 





140 


activity and an increase in the +H bay at COLL. BY 0821 UT 

(see Fig. 4.21B) the H perturbation at FITCH reached a value 

of -218y with the moment of the current system located 

between FTCH and FTSM. During the next two minutes the 

northern border of the electrojet rapidly developed resulting 

in a northward motion of the moment of the current System. 

Associated with this enhancement of the northern border, a 

large +D perturbation developed (as shown in Fig. 4.21C) 

indicating that a westward travelling surge was located over 

Our station line. By 0823 UT the maximum H perturbation along 

our line of stations reached a value of -374y and the moment 

of the current system was located slightly north of FTSM. 
Between 0823 UT and 0825 UT (see Fig. 4.21D) the 

H component at FTRL strengthened indicating an enhancement 

of the northern border of the electrojet. In addition, a 

negative H regime developed from FTMU southward as illustrated 

in Fig. 4.22. The +H regime that existed prior to this 

enhanced -H regime is also depicted in this contour diagram. 
The latitude profile for 0825 UT shows an example 

of the latitudinal asymmetries that occur in conjunction 


with the polar electrojet. It is interesting to note 
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the shape of the Z component south of the Z = O crossover. 
This deformation is similar to that shown in Chapter II for 
an asymmetric current distribution (see Fig. 2, Sas 

Immediately following the enhancement of the northern 
border, the pre-existing electrojet developed over the next 
few minutes forming a wide and fairly uniform current system. 
In Fig. 4.23 the latitude profile and polar plot for » 0831 UT 
are shown. The moment of the current system was located at 
~ 70.5° with a width of » 9°. The horizontal perturbation 
vectors shown here are indicative of a three-dimensional 
current system with a central meridian » 15° east of our 
line of stations. Note that the large perturbation vectors 
coincide with the Q = 2 oval. Also, the shape of the D 
profile is what we would expect to observe with a line of 
stations on the evening side of the oval (see Fig. 2.17 in 
Chapter II and also Chapter VI). 

Throughout the recovery phase of the substorm 
the current system gradually diminished in intensity with 
minor quasi-periodic enhancements at the northern border 
(similar to those observed during the substorm discussed 


in this section). 























eect? : ag 7 LA : 
. . ‘vb ae 

TS <r 7 a IVR, ved i 
Ve s 7 


a 
al 


-1evoz2019 0 = § ant to ‘doves saa ‘—on 


uo 


107 ou xgtqsnd nt nwone tent of xvelimte et not sr 
(48. S .pti 99a) not feurweses Ina vis tags 


nvadt yon. ant to snomegonsine ont pri twol [Tot vlosetbemal 
sxan sat avo bsqol gall batordodre. gisele 


madeve tn545B> mo thaw Yivrst bas sbtw 6 patriot eetu 


TU f€80 » tot Fofq +Sloq bas slrtorva sbusttsl ond es. re a i 
| oT 


7 


36 bstso0l 2sw meat2eye trex ww) Sdd to tnomom sdT -nwod 


oh 
hottedyvsveq [étnostyod AT .°R u to dtbiw 6 dttw ° OT 
a) ae a a 4 ‘a. 5 . a | 2 : ME 
lanorensmio-9s%it 6 to svitsotbnf ovs syed rietd eases 

he ~~ ie 

, 0 va 


‘- . 4 r r ‘ > a ee 

WO TO tess “el »v nerbitrem [stdnes 6 Atiw metaye, ne vw 
eyortasyv sorted TY OO an* lf ana+ ng gto an + ~ > . .% 
bw SOY FOSS UW. 90 SPp16i Sani tea3 830) -2no0rtsere ' anf 
— 5 i» 

G@ ott to eqed2e ott -o2elA .feve S = 0 oh¥ dttw eb onto 

, mm : 
§ ' , 


to eaff s dttw svysedo o+ tosqxs biuow ow tedw ef aft 


i © 
= 
Fon 


nr VI.S .phi 992) Tevo sAt to sbie pniinsvs oA? no eno seh: 


.(IV yvotqend o2fa Bas 1 xstase 
; mrosedue ont to sesdq \1svode, sit tuotpwor vt 


AS tw YStenstni nf bedetntmrb xi Teubsrp 0 moteve tnset 


i 





28 
"SO70d-Nishtion ont +5 2¢namsonsing ot Deron eaaes 


y os dune tb mrot edie ott. entvab ten anade s2ont OF 46 1 
-— i ( * 


—_ 










143 


S8 





"IN LEBO ~ sod (LPAO Z = D0) FOLd uelod pue aLLjoud apnziqey ezy “614 


SONLI Lb) 


08 s/s 








OL 


og 


O9 





‘Jas Oi 


“NIWOE YH 8B 


S61 ABO 





0 
ALISNSLNI 01314 


2 Ohl 
SEWN) 


08 


Och 


( 


“Jas ch 


"NINOE YH 8 SBI AYO 
SWIL IdSYSATNN 





06 


S3LYNIGYHOOD JIL3NObWO39 03193bH0 


O>z ---- 0 
O<z ---- + 


SYOLI3A NOTLYSYNLY3d JILINIWW TWLNOZIYOK 








_—_- ate 


~ 9. amit penvind 


7 A. 


alien MINCE AH 8 BI TAD 











. 
2) 





(CHE 





L 





* 
ba 
J 





EIErD IMLEv2 





~~ — i on 


rong ebuttiss 


ae S. | | 
= re 





144 


CHAPTER V THE POLAR MAGNETIC SUBSTORM AND THE ASSOCIATED 
AURORAL SUBSTORM 


Son. ANtroduction 


In this chapter we present the analysis of four 
polar magnetic substorms with the associated auroral sub- 
Storms. Since the auroral substorms were studied using data 
from only one all-sky camera, a full correlation between 
auroral and magnetic substorms was not possible. However, 
we were able to obtain very interesting results due to the 
fact that the all-sky camera at FTSM happened to be in an 
appropriate position to observe a major portion of the auroral 
activity during at least two of these substorms. During one 
of the substorms, that occurred on October 7, 1970, the all- 
sky camera was located west of the position where the west- 
ward travelling surge originated. Therefore we were able 
to study the magnetic effects associated with the passage of 
this surge over our line of stations. 

The substorm that occurred on September 1, 1970 
afforded us the opportunity to study in detail the latitudinal 
expansion of the westward electrojet within and probably close 
to the western edge of the auroral bulge. Also, during this 


particular substorm the position of the all-sky camera was 
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centered latitudinally with respect to the region of the major 
auroral activity associated with the recovery phase (see 
Akasofu, 1964). Recently the results of a study of the auroral 
substorm using a meridional line of stations in the Alaskan 
sector were reported by Snyder and Akasofu (1972) (also see 
Akasofu et al., 1971). It is quite obvious from the results 

of their study that analyzing an auroral substorm with data from 
several all-sky camera stations is necessary in order to 
minimize the possibility of producing erroneous interpretations. 
With this limitation in mind we present the following results 


of our study. 


5.2 Substorm events of October 7, 1970 (Day 280) 


Even though October 7, 1970 appears as an interna- 
tional quiet-day—COG)- with—-Kp+-indices-of—(0+,—-1- ,- 145; lo, lo, 
1-, O+, 1+; ZKp = 7-), a series of four substorms was observed 
in the Canadian sector between 0600 and 0900 UT. The H, D 
and Z components of the magnetic perturbations and Pi 2 micro- 
pulsations associated with three of these substorms are shown 
in Fig. 5.1 and 5.2 respectively. The symbols on the H compo- 
nent magnetogram allow the traces from each station to be 
followed, Prior togtene first™subdstorm shown in) Figs, 5.1, the 
development of a small eastward electrojet was observed flowing 


over our line of stations and centered between FTCH and FTSM. 
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At first we thought the development of this particu har acurrent 
System was indicative of a growth phase for the substorm that 
occurred at ~ 0659 UT. This was quickly ruled out when we 
discovered that a very isolated substorm (onset time of ~ 0615 
UT) was recorded at CHUR (see Ping. #5.3).) Pasitive devilatiions 
of the H component (with respect to the baselines ini Fig. 5.3) 
commencing at about the onset time of this isolated substorm 
are also discernable at ALER, RSLT, MOUL, BAKE, and BARR. 

The all-sky camera at FTSM did not show any auroral activity 
before 0659 UT. Also no Pi 2 activity was recorded with 


Our line of stations before 0659 UT (see Fig. 5.2). 


9.2.1 Polar substorm (0659 UT) 


The initial development of the first substorm 
presented in this section is illustrated in rig. Date amd ©. 
The plane projection of the all-sky camera data is included 
with each successive latitude profile. It should be noted 
that north (north being in the same direction as the H 
component of the perturbation field) is to the left in these 
diagrams. At 0659 UT a weak westward electrojet developed 
Slightly south of FTSM with a maximum H perturbation of 
« -20y (seesFiqg. 5.2 for the onset of Pi i2activity (/Asso- 


ciated with the onset of this weak electrojet was the appearance 
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of a bright E-W aligned arc located ~ 80 km south of FTSM. 
One minute later a westward travelling surge was recorded 
advancing from the east of our line of stations. The D pro- 
file became well deformed by the positive D perturbation 
associated with this advancing surge. The Shape of the D 
profile at 0659 UT is indicative of a current System with a 
central meridian east of our line of stations (possibly where 
the surge originated). At ~ 0701 UT the surge passed over- 
head moving to the northwest. Our prediction of the maximum 
H perturbation at this time is ~ -200y (note the change in 
the amplitude scale for each successive latitude profile). 
The surge continued in a northwesterly direction and by 0707 
UT passed out of range of the all-sky camera at FTSM. By 
0712 UT the moment of the current system was located at 

av 73.0°., the maximum northward motion that could be determined 
with our line of stations. The northern portion of the 
current system rapidly decayed thus resulting in a fairly 
rapid southward motion of the moment of the current system. 
During the recovery phase the moment moved to a latitude of 

ee dos (~ 2° north of the latitude at which the current system 
began) and by 0740 UT the substorm and hence the current system 
had almost completely disappeared. No more auroral activity 


was observed at FTSM after 0707 UT, therefore we could not 
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determine whether possible equatorward motions of auroral arcs 
were concurrent with the rapid southward motion of the current 
System during the recovery phase. 

The dynamic development of this particular substorm 
and the two later ones can best be illustrated by component 
contour diagrams (see Fig. 5.5a,b and c). The H and Z 
components from BAKE were used in the respective diagrams. 

The H component diagram shows the northward and southward 
motions associated with the expansive and recovery phases. 

We also see the +H regimes north and south of the current 
System as expected for a three-dimensional E-W current system. 
These are also evident in the latitude profiles shown pre- 
viously. In the region of CAMB the H component is slightly 
negative and may be due to the effects of polar cap currents 
as discussed in Chapter IV. Another interesting aspect of 
this diagram is the +H regime centered slightly south of 

FTCH occurring after the recovery phase. Associated with this 
positive regime is a +Z regime south of FTMU and -Z regime 
north of FTMU confirming the existence of an eastward elec- 
trojet. This is also apparent from latitude profiles and the 
component magnetograms (see Fig. 5.1). The buildup of an 
eastward electrojet following a substorm is not uncommon and 


may possibly explain the +H perturbations following -H bays 
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H component contour diagram for the series of substorms shown in Fig. 5.1. 
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Same as Fig. 5.5a but for the Z component. 
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reported by Heppner (1954). 

The D component diagram (Fig. 5.5b) is interesting 
in that it shows the northward motion of the +D regime 
associated with the surge propagating to the northwest. This 
is also apparent in the latitude PrOrvies alr iguwese4, The 
Characteristic latitudinal width of this and other similar 
+D regimes associated with surges and 100DS 15 4 400 km... If 
we assume a characteristic time of ~ 6 min from the contour 
diagram and the velocity of the surge to be ~ 1.1 km/sec (as 
obtained from all-sky camera data), then the longitudinal 
width is ~ 400 km. Current systems that possibly explain 
these +D component perturbations are presented in the next 
chapter. 

After the passage of the surge, the D component 
diagram shows the characteristics that one might expect to 
observe at a line of stations west of the central meridian. 
The enhanced -D regime at » 71.0° is believed to be caused 
by a NW-SE tilt of the electrojet. Polar cap currents appear 
to be the cause of the -D regime north of CAMB. In fact such 
polar cap currents usually affect all components in this region 
during most of the substorms we have analyzed thus far. 

The Z component diagram reveals most of the perti- 


nent information associated with the poleward expansion of 
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this current system. The average velocity for the northward 
motion of the moment of the current system for the first 9 
min was ~ 0.8 km/sec. During this period, the width of the 
electrojet increased rapidly from ~ 1 at 0659 UT to. 5 to 


re 


at 0710 UT. Although the southern edge (-Z peak) of the 
current system moved northward during the expansive phase the 
velocity was much less than that for the northern border. 

This would have to be the case since the width also developed 
rapidly during the same interval of time. In the next section 
we present an example of a poleward expansion with the sou- 
thern border remaining almost stationary. 

The development of the perturbation field as 
measured at 25 magnetic observatories is depicted by polar 
Blots shown in Fig. 5.6-and 5.7. Fig. 5.7 shows the hori- 
zontal magnetic perturbation vectors observed at » 0703 UT. 
The open circle to the east of our line of stations repre- 
sents the inferred location of the center (6.5) of the 
electrojet. This is easily verified by noting the distribu- 
tion of perturbation vectors south of the auroral (Q=1) oval 
(see Fig. 3.4). The large perturbation vector at FTSM is 
still strongly influenced by the perturbation field associated 


with the westward travelling surge. The surge probably 


originated south of the open circle. Fig. 5.7 shows the 
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distribution of perturbation vectors 5 min later (. 0708 UT). 
The D profile at this time had almost recovered from the 
distortion caused by the surge. Also at this time the -H 
perturbation reached its maximum value of ~ -450y as recorded 
with our line of stations. In both of these figures the | 
perturbation vector at BARR is probably representative of 
a high latitude eastward electrojet associated with the 
substorm current system. Also the rapid reduction in the 
amplitude of the Z component away from the current system 
in both latitude and longitude is very pronounced in these 
polar plots. 

Before the onset of the substorm there appears 
to be a slight depression of the H component at FRED, BOUL, 
DALL~ and TUCS {see' Figs-573). *This" could’ possibly’ be’ due 
to the asymmetric ring current associated with the eastward 


electrojet in the evening sector as discussed in Chapter IV. 


5.2.2 Polar substorm (0748 UT) 


The onset of the second substorm shown in Fig. 5.1 
was difficult to determine since it occurred west of our 
line of stations and the perturbation was of relatively low 
magnitude. Even though the onset of Pi 2 micropulsation 


activity was ~ 0748 UT (see Fig. 5.2), the first indication 
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of a westward electrojet was not observed by our line of 
stations until 0752 UT. The all-sky camera at FTSM recorded 
only 2 min of auroral activity (0753 and 0754 UT) during this 
entire substorm. At ~ 0753 UT a quiet arc was evident ~ 100 km 
south of FTSM oriented from WNW to ESE. The arc also coin- 
cided with the moment of the electrojet evident at that time. 
At 0754 UT only a diffuse auroral patch was observed ~ 150 km 
west and 50 km north of FTSM. 

The motion of the current system is clearly illus- 
trated in Fig. 5.5C. The velocity of the poleward motion 
between 0753 and 0759 UT was ~ 0.3 km/sec with the moment 
reaching 68.3° at 0759 UT. A very slow southward motion 
(probably associated with the recovery phase) started at 
~ 0759 UT and continued until the onset of the third substorm 
at v 0814 UT. 

Fig. 5.8 shows the distribution of the magnetic 
perturbation vectors observed in the northern hemisphere 
(see Fig. 5.3) along with the latitude profile for ~ 0800 UT. 
The D profile in this figure is an excellent example of a 
profile indicative of a current system with a central meri- 
dian west of our line of stations. This is clearly demonstrated 
io rigd.,.o.900a1so0%..Théescircle on-the-Q = |-oval represents 


the inferred center of the current system (» 20° west of our 
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Time of stations). ‘This could possibly explain the delay 
time (~ 4 min) between the onset of the Pi 2 micropulsations 
(~ 0648 UT) and the first Sharp perturbation of the electro- 
jet (0752 UT) as observed with our line of stations. 

The maximum width of the current system as observed 
by our Tine pf stations was «= 446-5.0° along with a maximum 
H perturbation of ~ -146y (0800 UT). However, these values 
could possibly be quite different from those that would 
have been observed if our line of stations had been located 


fartiner to the west. 


weews Polar substorm (0814 UT) 


During the recovery phase of the second substorm, 
the development of a third substorm started at a higher 
latitude (1 ToS") (see Fig. 5.1 and 5.5) and slightly west 
of our line of stations. From both the onset of Pi 2 micro- 
pulsation activity and the development of a westward electro- 
jet, the onset time was found to be »~ 0814 UT. The only 
auroral activity observed during this particular substorm 
was ate~v 0814 UT and 0815 UT... This took’ the form of a faint 
auroral surge located ~ 200 km north and 150 km west of FTSM 


moving rapidly to the northwest. 
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The development of this substorm is very similar 
to that of the previous one except for the fact that the 
electrojet was located at a higher latitude at the time of 
onset (x 70.5°). The velocity of the poleward motion of 
the Z = 0 crossover is about the same as for the previous 
substorm (~ 0.3 km/sec) with the moment of the electrojet 
reaching a latitude of %» 71.5° at 0820 UT. Commencing about 
the same time, a slow equatorward motion was observed for 
the duration of the substorm. 

A polar plot of the perturbation vectors and 
associated latitude profile for 0816:26 UT are shown in 
Fig. 5.9. The D profile may be indicative of a tilt of 
the electrojet as was noted in our discussion of the first 
substorm (the -D regime following the passage of the 
westward travelling surge). The Z profile was difficult 
to interpret owing to the lack of stations in the vicinity 
of the electrojet. This particular profile could have been 
drawn with more pronounced peaks. Although there is a 
definite +H response at lower latitudes associated with the 
onset of the -H regime to the north, the +H regime can 
partly be explained in terms of the recovery of the electro- 
jet from the second substorm due to the choice of baselines. 


However, we are not able to determine the relative magnitude 
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of each individual contribution. 
Three-dimensional plots of the substorms presen- 


ted in this section are presented in Appendix A6. 


5.3 Substorm event of September 1, 1970 (Day 244) 


On September 1, 1970 two large magnetic substorms 
were recorded with our line of stations. The H, D and Z 
component magnetograms and the associated Pi 2 micropulsa- 
tions are shown in Fig. 5.10 and 5.11 respectively. The Kp 
Por cessor this day ere, Git. 40. 50 At the Ao 3+: 

Kp = 27) and the Dst indices for the hours 5-12 consecutively 
are (3, 4, -2, -8, -19, -31, -27, -15) thus indicating an 
increase in the symmetric component of the ring current. 

Only an analysis of the first substorm is presented in this 
section since the center of the associated current system was 
only a few degrees east of our line of stations. 

From the Pi 2 micropulsations and latitude profile 
analysis the onset of the first substorm was ~ 0655 UT. The 
gradual development of an eastward electrojet preceded the 
onset of the substorm by ~ 1.5 hours-(see Fig. 5.10).From 
the available magnetic data no definite substorm activity was 
evident during this period prior to 0655 UT. Fig. 5.12 shows 


the approximate location of the eastward electrojet at ~ 0603 UT. 
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Magnetograms for the substorm that occurred 
On September 1, 1970. 
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5.11 Pi 2 micropulsation activity (40-150 sec) 


associated with the substorms shown in ee Peo ae ee 
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CORRECTED GEOMAGNETIC COOROINATES 


Fig. 5.12 Polar plot showing the distribution of pertur- 
bation vectors before the onset of the first 
Ssubstorm. The large arrowheads represent the 
probable location of the eastward electrojet. 
The dashed line represents C.G. midnight. 
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Note the lack of horizontal perturbation vectors in the post 
midnight sector (the dashed line represents C.G. local mid- 
night). This could possibly be due to the sparse distribu- 
tion of stations in this sector. In fact WHAL observed a 
Z perturbation of ~ 20y at this particular time. At the 
present time we do not know whether or not the development 
of this particular eastward electrojet in the evening sector 
represents the magnetic signature of a growth phase. 

Similar to the substorm event presented in the 
previous section, the H component is depressed at BOUL, 
DALL and TUCS just prior to the substorm onset (see Fig. 5.13). 
The diurnal component at HOLU and GUAM is too large to infer 


such a depression. 


5.3.1 Expansive phase 


In this section a detailed analysis of the develop- 
ment of the current system associated with the expansive 
phase of the auroral substorm is presented. First, the over- 
all dynamic development of the current system is discussed 
and secondly the very rapid intensifications of the northern 
border that produced the poleward expansion are examined. 

The expansive phase of the auroral substorm first 


became evident at 0656 UT when a bright arc appeared »~ 400- 
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500 km south of FTSM. Since the all-sky photographs from 
MEAN were obscured by clouds the auroral activity in the 
immediate region of the onset of the westward electrojet 

(% 64.3°) could not be determined. The ensuing northward 
motion and development of the system of arcs was very rapid 
and at ~ 0706 UT the northern edge of the bulge passed 
overhead FTSM reaching a latitude of » 70° Dy 0710 Qiy- 
Although it is difficult to study the poleward expansion of 
the auroral bulge with one all-sky camera, from the analysis 
of magnetic perturbations associated with this expansion 

we believe that 70° represents the maximum latitude reached 
by the poleward boundary of the bulge. 

To illustrate the rapid northward expansion of the 
westward electrojet during the first 15 min following the 
onset of the substorm, a series of latitude profiles are 
Shown in Fig. 5.14. The dashed line in each column of pro- 
files represents the location (64.3°) of the moment of the 
current system at the beginning of the expansive phase. The 
rapid northward expansion can easily be inferred by comparing 
the Z = 0 crossover with the dashed line in each successive 
profile. The D profiles immediately following the onset 
were very irregular and differed from that which is expected 


for an E-W orientated three-dimensional current system. We 
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believe that the deformations of these D profiles were 
caused by irregular auroral forms associated with the 
expansive phase of the auroral substorm. By 0702 UT the 
D profile began to resemble that expected for a current 
System with its central meridian to the east of our line 
of stations. 

Fig. 5.15 shows the Z component diagram for the 
first 16 min following the onset of the substorm. Three 
important aspects of the poleward expansion of the electro- 
jet are illustrated by this contour diagram: 

(1) The moment of the current system moved rapidly 
poleward. The velocity of this motion varied from ~ 280 m/sec 
(0656-0702 UT) to ~ 1.6 km/sec (0702-0704 UT). The southward 
motion of the moment began at about 0708 UT. 

(2) The southern boundary of the electrojet 
(represented by the dashed line through the maximum negative 
Z regime) remained very stable throughout the expansive phase, 
whereas the northern boundary moved poleward with a nearly 
uniform velocity of ~ 1.1 km/sec. This is also apparent in 
FIG. abe 14. 

(3) The width of the electrojet rapidly expanded 
in a very short period of time reaching a value of ~ 900- 


1000 km by 0709 UT. 
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Even though the all-sky photographs associated 
with the expansive phase were very difficult to interpret, 
we were able to extract what we think represents the impor- 
tant morphological features associated with the poleward 
motion of the auroral bulge. Fig. 5.16 and 5.17 show some 
examples of the interpreted arc structures along with the 
pertinent latitude profiles. Fig. 5.16 shows an intensified 
arc ~ 200 km south of FTSM along with the subsequent forma- 
tion of a westward travelling loop. At ~ 0702 UT a diffuse 
arc developed north of the arc containing the westward 
Propagating loop. One minute later (% 0703 UT) the whole 
region south of FTSM was overexposed due to an activation 
of the existing auroral arcs. By 0705 UT a well formed arc 
was observed overhead FTSM moving to the northwest. At 0707 UT 
this particular arc (now north of FTSM) intensified while 
continuing to move to the northwest (see Fig. 5.17). 

The tilt of the auroral arcs south of FTSM is 
interesting in that it may explain the large -D regime (due 
to the orientation with respect to our line of stations of the 
ionospheric current flow within the arc) also located slightly 
south of FTSM. 

To study the perturbations of the current system 
during the expansive phase, we used 10 sec resolution 


latitude profiles. Such an analysis revealed a very inter- 
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esting effect taking place during the northward expansion. 


fig. 5.10A Snows the initia development of the current 
System ~ 1 min after the onset. The current system at this 
time is probably less than 1.5° in width with a moment 
located at ~ 64.5° (note the change in the amplitude scale). 
This particular current system intensified over the next 2 
min, at which time the maximum H perturbation observed along 
our line of stations was ~ -240y. The width increased to 


ae oo 


with the moment moving northward to just slightly 
south of FTMU. The current system remained relatively 
unchanged during the next 2 min and at ~ 0701 UT the 

northern border of the preexisting current system intensified 
in conjunction with the activation of the auroral arc shown 
in Fig. 5.16. Fig. 5.18B shows the differential profile 
Showing the location of the current system associated with 
this intensification. The dashed line in each graph repre- 
Sents the location (64.3°) of the current system at the time 
of the onset of the substorm. The +H regime to the south 
represents an eastward electrojet (or a reduction of the 
westward flowing currents in this particular regime). The 
shape of the Z profile south of 64° rules out the possibility 
of explaining this +H regime totally by field aligned currents 


associated with an E-W three-dimensional current system (see 


a 
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Chapter II, Fig. 2.4). This particular deformation of the 
Z profile can be caused only by an ionospheric current. 
The shape of the D profile is what we would expect for a 
Current system with a central meridian east of our line of 
stations. 

A second intensification of the northern border 
occurred at ~ 0703 UT (this intensification was probably 
associated with the intense activation of the auroral arc 
overhead FTSM at 0703 UT. The resulting differential pro- 
file depicting this enhancement is shown in PU ee oa, 
Except for the slight deformation of the D profile at FTSM 
(probably due to a NW-SE orientation of the electrojet), 
this profile can easily be explained in terms of a three- 
dimensional current system. As a result of these successive 
enhancements of the northern border of pre-existing electro- 
jets, the overall system expands rapidly northward as 
previdis ly sshown in Fig. 5.14, 5,15 and also in Fig.. 5.17D. 

It is interesting to note that from the polar plot 
(for ~ 0708 UT) shown in Fig. 3.3 (Chapter III), the central 
meridian of the current system was calculated to be » 10-152 
east of our line of stations. The central meridians of the 
perturbation current systems associated with the enhancements 
of the northern border were also to the east of our line of 


stations. 






















. Al owas ~ % a 
‘ * : ie = #Y 


—? Pe aA 
; ; : | 9; PA ie i i 
: Y oO i ! iy Ror 
ont To nortsmveteb rsluotiied ar (d. ‘S. aif 4s 
,tasviwo > bisdqzonot ns x vino bezuss cn 52 91 7 
6 10t tosaxe biyow ow tsdw 27 afitord G ont, 5 0 oe 


to enti wo to ses9 nstbitem Tevdns> 6 Atty mote 


iad | _ 


| mie 


yebyod nysatyon sdi to nortsortte 2naint bnoss2 ae , 
Vfdsdovq 25Ww Aotteofttensint afAs) TU EQVO w +16) w9¢ 
“6 [svo%us sit to noftavitos natant ana Ahi 0228 

-o1g fetinaystttbh patsfuesy sAT .TU EOX0 +6 M273 beosnan 


ty 


.J8T.¢ .pri nt mwote ef tasmaonsdas etay antsatgab, mn 


MeT? $6 s/itowq 0 sit to nottsmyoteb Jniprie odd 193, $4 92x 


bar Aw « . +4 . ; : lye ' ; . Pe ay 
e(tetortosls sit to nottstnstyo 322 s oF a 79, 
fou 


-s971dy 5 to amyst nr benrslqxe $d yitess. we3 ef tt 


Pott i b x 
as } - on I « We Me > ~~ > ) «fj es 4 tT - ; a ’ eek Y 
SVieessoque saat: oO tives) 5 2A P eye tne two) [6 soote i Eta 
i + ne 
[ rte } ) = 


) fT 
“OFS 9919 ONTlerxs-siq to Wsby09d iamiisen ait to 2tnemsonsiins 


' ¢ P 
“fe i) ae 
26 biswiiyon \lbigsy ebnsgxe mateye [Tsyevo edd, 2tat 


Bess 

aw - . * r ‘ r - . r . JP mits ns 

-G\i.c .Pfi of o2is5 base @l.e .B{.2 .ptt nt awode ylevotyes 
avin OF 


toig welog sit mov? tant ston ot pnttzevstnt et 41 fos 









7 : Pigs Ap. 2 
: fevtngo sat ,.(1iI yaegs 9) .£ £ .pi. nt nwode (TU 
save 


sp) Se, 


W2 tneriw2 s 






efor x od ot bate tuotes 26W i 


a 
0 enstbt bys Tevspe3: 


'7 
a be 


y & S37 BTOC eT: 2 
sake pozas amos: 


+ ey beg abe te 
Va Melina @ 


Ve 


: i 1 
Ne, tne" sue 
- te - a = 


183 


5.3.2 Recovery phase 


The beginning of the recovery phase for this parti- 
cular substorm is not as clearly defined as that for the 
0659 UT substorm (Day 280) discussed in the previous section. 
Even though the first equatorward motion of auroral arcs was 
observed at ~ 0710 UT, Snyder and Akasofu (1972) reported 
that such motion can occur while the auroral bUiae -18.561.11 
expanding. We also noted previously (see Pidse 5. love that 
the moment of the electrojet began to move southward begin- 
ning at ~ 0708 UT. To make the problem of choosing an onset 
time for the recovery phase even more difficult, the northern 
border of the electrojet experienced two major tices eres 
at ~ 0717 UT and at 0806 UT. In fact, quasi-periodic intensi- 
fications of the northern border were in progress throughout 
the entire substorm. Snyder and Akasofu (1972) also questioned 
the meaning of the recovery phase in terms of whether it 
represented the reappearance of a significant equatorward 
component of motion within the bulge, or the recovery of the 
poleward boundary of the bulge. For the purpose of the 
following discussion we set the beginning of the recovery 
phase at the time of the first southward motion of the 


moment of the electrojet (~ 0708 UT). 
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The overall dynamic development of the H component 
of this particular substorm is shown in Fig. 5.19 (3-D 
diagrams for the D and Z components are included in Appendix 
A6). A ‘double' current System, represented by the two 
‘mountain ranges' (centered about » 64° and 68°) in the 
3-D diagram, existed throughout the entire substorm. In 
fact a triple peaked current distribution was evident at 
~ 0717 and 0809 UT following northern border enhancements. 
The northern border enhancement at 0717 UT can not be seen 
in this particular perspective view of the H component array 
due to the strength of the second maximum at » 68°. 

The two polar plots shown in Fig. 5.20 show the 
distribution of the magnetic perturbation vectors at ~ 0730 UT 
and 0809 UT during the recovery phase. It is interesting to 
note the similarities between these two distributions even 
though they are separated in time by ~ 40 min. The longi- 
tudinal expansion of the electrojet can be seen by comparing 
these polar plots with that shown in Chapter III, Fig. 3.3 
for ~ 0708 UT. 

During the recovery phase, most of the auroral 
activity and perturbations of the current system were centered 
about FTSM. In fact, quasi-periodic intensifications of the 


northern border of the current system were in progress much 
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of the time. It was found that the D component responded 
to all the eastward and westward propagating auroral forms. 
These D component responses were similar to the +D perturba- 
tions discussed in the previous section for the westward 
travelling surge. To illustrate these responses, a filtered 
version (60-1000 sec) of the D component magnetogram is 
shown in Fig. 5.21. The dotted and solid lines represent 
westward and eastward propagating auroral forms respectively, 
that could be easily identified in the all-sky camera photo- 
graphs. During the first hour of the recovery phase these 
+D perturbations occurred at ~ 7 min intervals. The lines 
were drawn along the peaks of the +D perturbation maxima (in 
latitude) to emphasize the fact that such perturbation fields 
may have their lines of N-S maxima slightly tilted with 
respect to our line of stations. 

Upon closer examination of the Pi 2 micropulsations 
for this substorm (see Fig. 5.12), we noticed that an enhance- 
ment of Pi 2 activity was also associated with these eastward 


and westward motions of auroral forms. An example of a 
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perturbation attributed to the passage of a westward travelling 
Toop is shown in Fig. 5.22. The top profile (~ 0716 UT) 
represents the unperturbed system whereas the bottom profile 
(~ 1 min later) shows the perturbed system in response to 

the passage of a westward travelling loop. Note the deforma- 
tion of the D and Z profiles and the enhancement of the H 
component at the northern edge of the current system. When 
we set the baseline at the time of the first profile (0715:51 
UT), the resulting differential profiles shown in Fig. 5.23 
were obtained. The most apparent property of these differen- 
tial profiles is the +D regime similar to that normally 
observed during the passage of a westward travelling 

Surge. The second property is the occurrence of an apparent 
shear system of currents involving a westward electrojet 
north of the auroral form and an eastward electrojet to the 
Soutn (see Chapter 11, Fig. 2.18C). 

Fig. 5.24 shows the differential profile represen- 
ting the magnetic perturbation associated with an eastward 
travelling auroral form (showing the accompanying auroral 
diagram). In this particular example, the +D perturbation 
is not as large as those shown in Fig. 5.23 due to the fact 
that the main N-S aligned portion of the auroral form was 


to the east of our line of stations. At ~ 0833 UT, however, 
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the profile (not shown here) showed a much more enhanced D 
perturbation. 

We analyzed several cases of E-W motions of isolated 
auroral forms and found similar magnetic perturbations along 
our line of stations. We also discovered that this form of 
differential profile was associated with a slowly moving 
N-S aligned auroral arc which occurred on September 19, 1970. 
A model current system that may explain these perturbations 
is discussed in detail in the next chapter. 

The last major enhancement of the northern border 
commenced at ~ 0806 UT and the associated latitude profile 
(v 0810 UT) is shown in Fig. 2.25. At 0800 and 0801 UT a 
well defined westward travelling loop was observed ~ 150- 

200 km north of FTSM (the +D response is clearly evident in 


Fig. 5.21). Approximately 5 min dater a diffuse arc, located 


~ 300 km north of FTSM developed (0806 UT). During the 
following 3 min a rapidly propagating auroral form was 
observed to develop while moving westward. It is not certain 
whether this represented a westward travelling surge (as we 
have indicated in Fig. 5.19) or a propagating loop. A west- 
ward propagating auroral form was also observed during the 


enhancement of the northern border at » 0717 UT. 
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At ~ 0810 UT the resulting current system was 
~ 10-11° wide with both the Z and H component profiles 
shown in Fig. 5.25 indicative of complex structure in the 
latitudinal current distribution (possibly a triple peaked 
distribution). 

Throughout the remaining portion of the recovery 
phase the northern border decayed while the southern border 
and the portion of the electrojet around FTCH (67.0°) 
experienced marked fluctuations. By ~ 0930 UT (the time 
of the onset of the second substorm), the maximum pertur- 
bation observed by our line of stations was approximately 


150y. 
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CHAPTER VI MODELLING OF THE CURRENT SYSTEM ASSOCIATED 
WITH THE POLAR MAGNETIC SUBSTORM 


SE LRLroduction 


In Chapter II and Appendices Al, A2 and A3 we 
presented the results of our study of magnetic field per- 
turbations generated by three-dimensional current systems 
along with a convenient method of evaluating these fields. 

In this chapter a comparison is made between 
the observed magnetic perturbations associated with polar 
magnetic substorms and the magnetic fields associated with 
three-dimensional current systems. 

In section 6.2 we examine the best possible 
depth at which to place the superconductive layer below 
the surface of the Earth in order to compensate for 
PoUeT i iOng Bi fects. .Section.6.3. contains a discussion, of 
coordinate systems along with the presentation of the 
results of our modelling of the substorm associated 
current system. The modelling of magnetic field pertur- 
bations associated with eastward and westward propagating 


auroral forms is discussed in Section 6.4. 
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6.2 Depth of the superconductive layer used in our modelling 


Since there are many parameters which affect the 
Shape of the latitude profiles associated with model three- 
dimensional current systems, we decided to incorporate the 
results of previous research in order to remove one of these 
parameters, namely the depth of the superconductive layer. 

McNish (1938) was the first researcher to examine 
the effects of induction on the magnetic field perturbations 
due to polar electrojets. In his model he ignored the 
curvature of the Earth and represented the electrojet by 
an infinitely long sheet current of finite latitudinal 
width, positioned at a given height above the surface of 
the Earth. To allow for the effect of induction he included 
a superconductive layer at some depth below the surface of 
the Earth. With this particular model McNish was able to 
obtain a fairly good approximation of the observed data for 
a particular event by assuming a sheet current of » 670 
km wide positioned at a height of 100 km, and a super- 
conductive layer 200 km below the Earth's surface. Forbush 
and Casaverde (1961), using a similar model, studied the 
magnetic field perturbations due to the equatorial electro- 
jet. They found that a depth of 250 km for the position 


of the superconductive layer yielded the best re- 
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sults. They also applied the method of separating the 
external and internal fields used in their study of the 
equatorial electrojet, to the auroral zone magnetic data 
analyzed by McNish (1938). The depth they determined for 
the top of the superconducting layer was 250 km, in complete 
agreement with that obtained from the study of the equatorial 
electrojet. Czechowsky (1971), using the same model but 
allowing for various latitudinal distributions of current, 
obtained an average value for the depth of the supercon- 
ductive layer of ~ 450 km. In his study, however, Czechowsky 
used magnetic data from only four stations, which as we have 
shown, may produce ambiguous results in terms of the causative 
current system. 

Since the equatorial electrojet is probably more 
representative of a two-dimensional current system and also 
a more static system, we decided to use the value of 250 km 
for the depth of the superconductive layer as obtained by 


Forbush and Casaverde (1961). 


6.3 Modelling of the polar electrojet 


In Chapters IV and V we used the corrected 
geomagnetic (C.G.) coordinate system in our analysis of 


the dynamic development of polar magnetic substorms. Since 
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this particular coordinate system is non-orthogonal when 
transformed into geographic coordinates, distances measured 

in degrees of C.G. latitude and C.G. longitude are not 

the same as those in geographic or centered dipole (geo- 
magnetic) coordinates. Therefore, since the model current 
systems are based on an orthogonal system of coordinates, 

we used the geomagnetic coordinate system for the presen- 
tation of the observed magnetic data. To do so the hori- 
zontal components of the magnetic perturbations were trans- 
formed from HDZ coordinates to H'D'Z (geomagnetic) coordinates 


using the following equation 


H' cos 8 -sin 0 H 
= bis 1 


D' sin 9 cos © D 


where 6 is defined as the difference between the magnetic 
declination (¢) and the dipole declination (w) as 
illustrated in Fig. 6.1 (X and H' being in the direction of 
geographic and dipole (geomagnetic) north respectively). 
Table 6.1 contains the values of 90, 6 and w along with 
the geomagnetic coordinates for each observatory whose data 


were used in this thesis. 
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DIPOLE 





Fig. 6.1 The relationship between the coordinate systems 
defined within the text. 


Three of the latitude profiles we have modelled 
thus far are presented Tn Fig. 6.2, 6.3 and 6.4. The first 
two profiles are for the substorm that occurred on June 15, 
1970 (Day 166) and the remaining one is for the substorm that 
occurred on September 1, 1970 (Day 244). The parameters 
for the three model current systems we believe to be re- 


sponsible for the observed magnetic perturbations are listed 
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INEORETLICAL PROFILE 


Me 5 oh ul MIN. 2ecet . 
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Comparison of theoretical and observed profiles for 
Day 166 at ~ 0715 UT. The arrows represent the 
boundary of the oval along the central meridian. 
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Pig. 6.3 Same aS Fig. 6.2 but for Day 166°at ~~ Oy 43eUT- 
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Fig. 6.4 Same as Fig. 6.2 but for Day 244 at » 0704 UT. 
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in Table 6.2. The oval parameters were discussed previously 
in Chapter II (see eq. 2.25 and Fig. 2.15). The model 
systems for 0715 UT and 0743 UT show the expansion of the 
Current system associated with the development of the 
northern border. During the interval of time represented 

by these profiles, the current system expanded both in 
length (32.5° to 50°) and in width (3.5° to 6.5°). The 
total integrated current increased by a factor of 3, from 
0.23 x 10° amp to 0.69 x 10° amp. 

In all three profiles, we see the deformation of 
the D' component due to the fact that the ionospheric portion 
of the current system was flowing along the auroral oval 
rather than along a constant latitude circle. Also evident 
in these figures is the relatively good fit of the 
theoretical Z profiles to the observed Z perturbations. 

It appears that the choice of 250 km for the depth at 
which to put the superconductive layer was satisfactory. 

The polar plots for both the observed and 
theoretical magnetic perturbations for these times are 
Shown in Fig. 6.5, 6.6 and 6.7. The expansion of the 
current system for Day 166 is easily seen by comparing 
Fig. 6.5 and 6.6. The eastward electrojet in the evening 


sector and the currents in the polar cap region probably 
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account for the disagreement between the observed and 
theoretical perturbation vectors in these regions. The 
polar plots shown in Fig. 6.7 are remarkably similar, 

even in the polar cap region (except ALER which is extremely 
anomalous (see Whitham and Andersen, 1962)). The substantial 
disagreement for some stations close to the electrojet can 
easily be reconciled by the fact that latitudinal gradients 
of the H component in this particular regime are extremely 


large. 


6.4 Model current system to explain the magnetic perturba- 


tions associated with E-W propagating auroral forms 


In Chapter V we presented examples of differential 
profiles showing the magnetic perturbations associated with 
eastward and westward propagating auroral forms. In this 
section we present a current system that could possibly 
explain the observed perturbation patterns. 

A series of three-dimensional current configura- 
tions are presented in Fig. 6.8 to illustrate the changes 
that likely take place during the passage of an auroral form. 
Configuration A represents the unperturbed system and D the 
perturbed system. The southward directed current (i,) 


indicates the location of the N-S aligned portion of the 
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auroral form. To model this particular system, we introduce 

a shear of A to produce B. We then add B to the narrow 

N-S three-dimensional current system C to produce D. We 

are not suggesting, however, that this sequence is necessarily 
representative of the dynamic development of the configuration 
D. The differential profile (D-A) as observed along the 
meridian ¢ due to the change from configuration A to D, 

is shown in Fig. 6.9. The differential profile (B-A) is 

also included in this figure. The original system (A) is 


a three-dimensional E-W current system with the following 


parameters: 
I a ese ieee 2 ao pe 40° 
Cea? re Serer 5° 
a) 
86 | ere Sere ee ee ee 61.45 
Total ‘integrated currenti:.......0.- I 0 10° amp 
Depth of superconductive 
ee dad) a SP eee were et ee ee ee 250 km 


The shear is along the central meridian with a i displace- 
ment; thus the N-S current (C) is 0.2 ~x 10° amp. 

Although we have not attempted to model a specific 
event, it is interesting to note the similarities that exist 


between profile (D-A) and the example SHOWN AT iF POs Secon 
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This particular model could also be used to explain the 
magnetic perturbations due to the westward travelling surge. 
For example the current configuration D may develop along 
with an enhancement of the northern border of the eastern 
portion of the sheared system. This could be associated with 
outward flowing Birkeland currents at the northern border near 
the shear zone (o0) and an enhancement of inward flowing 
Currents near the meridian ($-). 

Other possible current configurations which should 
be considered when modelling perturbations of current systems 
in the polar region are shown in mg. «6b. ion fOY a N=5 
three-dimensional current system 60° in length and 4° in 
width with a total integrated current of 2.0 x 10° amp, the 
observed perturbation along the central meridian is shown 
in Fig. 6.11 by the D profile in each latitude profile. 

The H and Z components of the field are zero along this 
meridian. If we introduce a shear along the central meridian, 
the H and Z components will no longer be zero. The latitude 
profiles for the four different current configurations in 

Fig. 6.10 are shown respectively in Fig. 6.11 (no induction 
field is included in these profiles). From these profiles 

we see that apparent enhancements of both eastward and west- 
ward electrojets could possibly be due to shears of N-S 


three-dimensional current systems. 
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CHAPTER VII DISCUSSION AND CONCLUSIONS 


7.1 Morphology of polar magnetic substorms 


In Chapters IV and V we presented a detailed study 
of the dynamic development of polar magnetic substorms. Using 
the framework proposed by Akasofu (1964) for the development 
of the auroral substorm, we now present what we feel to be 
the important morphological features of the polar magnetic 
substorm as follows: 
(i) The intensification of the westward electrojet 
near the southern border of the auroral oval is coincident 
with the brightening of the southernmost auroral arcs 
(expansive phase Stage I). 
(ii) The northward expansion of the northern border 
of the electrojet may take place in either a continuous or 
‘step-like’ manner. During this poleward expansion the southern 
border of the electrojet usually remains relatively stationary 
although it may advance poleward with a velocity which is much 
less than that of the northern border. This development of 
the electrojet is associated with the expansive phase Stage II. 
(iii) An intense enhancement of the current flow at the 
northern border of the electrojet along with the development of 


a +D component regime which is confined in both latitudinal and 
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longitudinal extent. This occurs in conjunction with the 
westward travelling surge during expansive phase Stage III. 
The +D component regime is believed to be caused by a north- 
to-south current flow in the region of the surge. 

(iv) Associated with the recovery phase Stage I, 
minor quasi-periodic intensifications of the northern border 
of the electrojet occur in conjunction with eastward and 
westward propagating auroral forms such as loops and 
irregular folds. These intensifications can also occur at 
other latitudinal positions and are responsible for the 
variations in electrojet intensity pointed’ out: byc Kamide 
et al. (1969). 

(v) Equatorward motion of the moment of the 
electrojet due for the most part to the decay of the 
Current system at the northern border. This is associated 
with the equatorward drift of the auroral arcs during the 
recovery phase (Stages II and III). 

The observed pattern of substorm development 
discussed above leads us to the following suggestions 
regarding the morphology of the polar magnetic substorm. 

We first note that the poleward motion of the 
electrojet may occur in "discontinuous" steps. 


Such a pattern of ‘development has» recently ° 
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been suggested by Hones et al. (1972). That is, rather 

than having the main electrojet expand northward, new current 
elements may suddenly appear to the north of the moment of 
the main electrojet. Thus the description of poleward 
expansion by Akasofu et al. (1966b) must be treated with 
discretion in that it implies that the motion is continuous. 

We have found that the first major enhancement of 
the northern border of the electrojet is associated with the 
generation of the westward travelling surge. The large 
positive D perturbation associated with this intensification 
is easily identifiable on normal magnetograms, and hence 
can be used to identify the presence of westward travelling 
Surges. 

The development of the westward surge (expansive 
phase Stage III) marks the onset of an interval of quasi- 
periodic intensifications of the northern border of the 
electrojet (recovery phase Stage I). The morphology of all 
these northern border disturbances seems very similar, so 
that we feel that expansive phase Stage III and recovery 
phase Stage I are part of the same process and should be 
combined under a common phase. We have suggested (Kisabeth 
and Rostoker, 1971) that this phase involves the direct 


input of energy from the magnetotail in impulsive bursts. 
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Finally we note that the main electrojet is often 
very stable during the period of northern border activity 
following the initial poleward expansion. Therefore we 
suggest that this electrojet is driven by a different source 
mechanism than the northern border activity. 

It should be noted that this pattern of develop- 
ment of the westward electrojet is supported by the theo- 
retical studies of Coroniti and Kennel (1972) who claimed 
that the onset of substorm activity was initiated through 
resistive instabilities at the inner edge of the plasma 
Sheet which maps to the equatorward edge of the oval. 

The development of magnetospheric substorm activity 
appears to be associated with the growth of an eastward 
electrojet in the post-noon and evening sector. This eastward 
electrojet can be very broad and stable, even during periods 
of substorm activity (see Rostoker, 1972). However, we have 
shown in Chapter IV that marked intensifications of the 
eastward electrojet may occur in conjunction with the onset 
of substorm activity. Therefore we feel that the eastward 


electrojet does respond to substorm fluctuations, although 
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it is not presently possible to evaluate quantitatively the 
level of response. We do feel, however, that this indicates 
that a single three-dimensional current system such as that 
Proposed by Akasofu and Meng (1969) and Bonnevier et al. 
(1970) is insufficient to describe the overall substorm 
perturbation pattern. It is probable that the combination 
current system described in Chapter IV (see Fig. 4.1) 
coupled with a realistic polar cap current system (or 
equivalent field aligned current flow) will eventually be 
necessary to completely describe the observed perturbation 
pattern. 

Finally we wish to point out that this study has 
emphasized how difficult it is to define the onset time of 
substorms from normal magnetograms only. In many cases, 

a strong perturbation associated with a northern border 
intensification in the middle of a substorm event, appeared 
as an "isolated" substorm onset at auroral zone stations 

far from the intensified portion of the auroral electrojet. 
The localized character of many of the substorm perturbations 
confirms the observation by Rostoker and Camidge (1972) that 
azimuthally confined sectors of the magnetotail may be 
sequentially activated thus generating the quasi-periodic 


activity at’ the northern border of the electrojet. 
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7.2 The modelling of polar magnetic substorm 


Current systems 


In modelling the polar electrojet, we have found 
it necessary to develop much more sophisticated models 
lay have been presented by other workers (eg. Akasofu 
and Meng, 1969; Bonnevier é@t al., 1970; Kamide and 
Fukushima, 1971). In particular we have included the 
following variable parameters in our model current systems: 
(i) A finite length and width of the ionospheric 
part of the current system. 
(ii) Asymmetric flow of Birkeland currents along 
the length of the electrojet. 
(iii) Combination current systems involving both 
eastward and westward electrojets. 
(iv) Various latitudinal current distributions. 
(v) Flow of current along the auroral oval. 
(vi) The effects of induced currents in the Earth. 
We have found that the variation of these 
parameters can significantly alter the computed perturbation 
pattern along a geomagnetic meridian for the three-dimensional 
system. To provide more information as to what these 
parameters should be, we used supporting data from the world 


wide network of observatories as well as all-sky camera data. 
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The examples of the modelling of substorm associated 
current systems given in Chapter VI clearly demonstrate the 
viability of our modelling techniques. It should be 
emphasized that the three-dimensional models we use repre- 
Sent only ‘equivalent’ current systems and should be treated 
aS SUCTT: 

Our study of induction effects substantiates the 
findings of Bonnevier et al. (1970) who concluded that the 
practice of correcting the observed components by constant 
factors (2/3 for H and 2 for Z) can not be used in studying 
localized currents such as the polar electrojets. Therefore, 
the importance of including the effects of induction in any 
model studies cannot be emphasized enough. For example 
Kamide (1970) presented variations of a model Current system 
which included ionospheric return currents rather than field- 
aligned currents in order to explain the observations of 
Walker (1964) and Bonnevier et al. (1970). If the effects 
of induction had been included in this model system, the 
relative strengths of the H and Z components would have been 
unrealistic. 

Finally, until now the expense of evaluating the 
magnetic effects associated with complex current systems 


has been insurmountable. - Since’ the>techniques: devetoped in 
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this thesis allow such evaluations to be done with relatively 
small amounts of computer time, we feel that these techniques 


can be incorporated into any routine data analysis. 


7.3 Techniques of presentation of polar magnetic 


substorm data 


Throughout the course of this study, the various 
techniques of data analysis and presentation were utilized 
(see Chapter III, section Su Jie ol sthis.section we present 
a brief summary of the attributes of each method. 

(i) Latitude profile: This method is very useful 
in providing a comprehensive picture of the distribution 
of the current in the westward electrojet at any instant 
during a substorm. A cursory examination of a profile can 
give considerable~ information about the position of the 
station line relative to the ionospheric portion of the 
current system, as well as a good quantitative estimate of 
the strength of the substorm. However all information about 
the temporal development of the substorm, as well as the 
longitudinal variation of the current system, is Tost in 


this presentation. 
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(ii) Polar Plots: This commonly used technique was 
modified to incorporate some Z component information and the 
auroral oval contour as a reference frame in which to view 
the data. Polar plots are particularly useful to gain a semi- 
quantitative idea of the world wide character of the substorm 
perturbation pattern at a given instant in time, along with 
Some indication of the azimuthal extent of the intensified 
portion of the electrojet. Quantitative information regarding 
the Z component perturbation is lost; however, our modification 
allowed changes in Z component polarity to be traced, yielding 
information on the position of the moment of the current 
System. As in the case of the latitude profile, information 
regarding the temporal development of the substorm current 
system is lost. 

(i171) Component Contour Diagrams: This technique, 
developed by Zaitzev and Bostrom (1971), is very useful in 
tracing the temporal development of the substorm. Information 
is presented showing the motion of the moment of the current 
system, and its latitudinal expansion and contraction. 
Quantitative information regarding the strength of each 
component is readily available, and polarity changes are 
well portrayed. The information missing in this format is 


the instantaneous world wide character of the substorm current 
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System. The method is also cumbersome in that it requires 
three separate diagrams (one for each component) to transfer 
the information completely. 

(iv) Three-dimensional diagrams: This technique is 
excellent in its ability to portray, using the technique 
Of perspective viewing, the temporal development of a polar 
magnetic substorm. It is particularly useful in illustrating 
the changes in the latitudinal current distribution during 
the course of the substorm as well as effectively presenting 
qualitative information on the variations in the strength of 
the substorm current system. Once again, information is lost 
in the world wide character of the substorm current system, 
and the method is cumbersome in that it requires three 
separate diagrams to transfer the information completely. 

From our study, we have concluded that no one 
method of data presentation adequately describes the character 
of the development of the polar magnetic substorm current 
system. Future studies should involve a combination of the 
techniques discussed above in order to yield complete 


information on the character of the substorms under study. 
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APPENDIX Al METHODS OF EVALUATING THE MAGNETIC FIELD 


PERTURBATIONS OF MODEL CURRENT SYSTEMS 


Al.1 Method using Biot-Savart's law 


In this section we apply Biot-Savart's law in 
spherical coordinates to an arbitrary current loop c 
exvernai to the Earth's surface (See Figé Al.1).« The 
result of this application isin the form of an integral 
of a simple third order matrix which is easily applied to 
an arbitrary three-dimensional current system in the 
magnetosphere. 

The differential magnetic field vector dBy at 
the observation position (1 2852¢,) due to an elemental 


current vector Jds is given by Biot-Savart's law 


ds x R 


ities fee 
|r| 


(Al.1) 


where 


Ran =vertorer (Aae2) 


and k is the proportionality factor dependent upon the 


system of units used. The field vector dB can be 
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SURFACE 


PIO ALA Diagram defining vectors in spherical coordinates 
for magnetic field calculations. 
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written in component form as 


dB = dBir_ + dB,o + Bad, (A1.3) 


Since the coordinate directions (7,654) at the source 
position differ from the coordinate directions (r58, 945) 
at the observing position, ds and B must be expressed 
in terms of (1385545) for “Als3i to be valid. If we express 
both ds and R in terms of (r,6,4) then the right hand 


Side of equation (Al.1) must be transformed as follows: 


(A1.4) 


where A is the orthogonal transformation matrix defined 


> 
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ar t 

8 x A |6 (A1.5) 
O 

7 " 


In the present derivation we express ds and R in terms of 
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(r 28, 2¢)) rather than transforming the resultant cross- 
product as illustrated.jin Al.4. To accomplish this we must 


~ 


still make use of the transformation matrix A 

We can derive A by making a series of two 
coordinate transformations; a transformation from (yr ,8,0) 
coordinates to (x,y,z) coordinates followed by a transfor- 


mation from (x,y,z) coordinates to (r 28024) coordinates. 


The first transformation is given by 


X . 
y : crate (A1.6) 
z $ 
where 
Sin 8 cos cos 8 tos @ -sin o 
Sy = sin 0 sin 6 cos @ sin 6 cos Gabe 
cos 6 -S71N» 6 0 


and the second transformation by 
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oA x 
85 = So y (A1.8) 
6 z 
where 
sin 85 cos uP sin Fs Sin 25 cos oe 
So = cos 85 cos b> cos 85 sin b5 -sin u 
-sin b5 cos 6 0 
(A1.9) 
From Al.5, Al.6 and Al.8 we obtain 
r r 
0 
Hf - S 5 Sy 6 
Ys w) 
(A1.10) 
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Therefore A is given by 
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A = Ao) 

Jo 

where 

Ay, a ein 8 sin 8 cos (4 
Ay» = sin 8, cos-9 cus (6 
Ay3 = sin 6, sin (6, - 6) 
An) = cos 0. Sin 6 cos (6 
Ano = cos @, cos 6 cos (6 
An. = cos 6) sin (>, - ¢) 
Ax) = -sin 6 sin ($5 - >) 
Ax» = -cos g sin (6, - ¢) 
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Ax, e- cos. (6. =4) 
VALS 2) 
The vector R= can be expressed as 
R = Piaget ee uet os (A113) 


whereas the vector ds is generally expressed in terms of 


the source coordinate directions as 


ds = dsyr + ds 0 fe ds ¢ (A1.14) 
Therefore, before the cross-product (ds x R) can be 
evaluated, we must transform ds to the (7, 29,2%,) 
coordinate system by the equation 
r 
0 
dg = [ds ds acta enn le (A1.15) 
] 2 3 ) 
wy) 


where A is the transpose of A. Thus 
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+ (dsj A., + ds, Az, + ds,A..) 6 (A1.16) 


The components Ry “aS and R in Al.13 can be derived 


by making use of ql as follows: 
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0 
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merefore; substituting Al.16 and Al.18 inte Alva, dR can 


be written in matrix form 
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gp #4, Oty, - (A1.19) 
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where 
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-rAj ods. (rA33 ~ rAjo)ds4 -(r Ang + rAz5)ds. 
(A120) 


Using Al.3 and A1l.19, the total contribution to the j 
field component at the observation position (r 28) 9 ) due 
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Al.2 Method using infinitesimal] magnetic dipoles 

The specific problem of eValuating the scalar 
magnetic potential and thus the magnetic field due to 
three-dimensional current systems which include dipole 
field-aligned currents has been solved using infinitesmal 
magnetic dipoles (Bostrom, 1968; Bonnevier et al., 1970). 
In this section we give only a brief summary of their 
solution to the problem. 

Bonnevier et al. (1970) considered the two current 
loops shown in Fig. Al.2 where field-aligned currents in 
both loops have infinitesmally small separation at the 
base of the ionopshere. Hence the ionospheric portion of 
the E-W and N-S loops have lengths of r,sin 8,do, and 
rjde, respectively. ((r, 0,5) is the coordinate posi- 


tion of the endpoints of both loops in the ionosphere). 





Fig. Al.2 E-W and N-S current loops with associated 


infinitismal magnetic dipoles (after Bonnevier et al. onl S204 
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The current in each loop is defined as 


dJ = Le (8) 56,) r,d6, 


dJ = i (6159, ) r, sin 6, do, (A122) 


where I and I represent the eastward and southward 
component of the height-integrated ionospheric current 
density if respectively. They have shown that the scalar 
magnetic potential due to these loops can be represented 
by an integration over elementary magnetic dipole moments, 
d?m , distributed along field lines as shown in Fig. Al.2. 
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d*m is defined by the equation 
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The coordinate system used is the same as that shown in 
bade Al.sle The Superscripts refer to the order of the 


differential. 


Since the scalar magnetic potential at (r 26 so.) 
we 


due to moment d°m at (r,6,¢) is given by 


Ane armel mies (A1.24) 


this expression can be integrated with respect to 6 
(integration along field lines) to obtain the total poten- 


tial due to both loops 


t/2 
di“e = k oe (A1.25) 
hid 


The vector R is the same as the one used in the previous 
section but is defined in terms of unit vectors r,06,¢ 


“n~ “n~ 


rather than Pasta sot . The components of R defined by 
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R = Ryr + R.8 + Ro (A1.26) 
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can be derived using the same method as was used in Al.1. 
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The magnetic field is obtained by taking the 


gradient of the potential, hence, 
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Therefore, for any distribution of current in the 
ionosphere, equation Al.28 can be integrated with respect 


to the ionospheric variables 6 and >, to obtain the 


1 
total contribution to the magnetic field at any observating 
position CBee Gad Wy This three-dimensional integration 
represents a volume integration over a distribution of 
magnetic dipole moments d?m . It is interesting to note 
that a fourth integration can also be performed with respect 
to ry provided f is expressed in the more general form 


Tt (r,.8,56,) - ‘thus an arbitrary distribution of current 


with height can also be considered. 
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APPENDIX A2 EVALUATION OF THE MAGNETIC@FIELD DUE TO 


CURRENTS INDUCED IN AN INFINITELY CONDUCTING 
EARTH 


We. | Introduction 


The derivation presented in the next section of 
this appendix is based on the results obtained by Ashour 
(1971) who calculated the magnetic field produced by currents 
induced in a uniformly finite conducting spherical shell 
where the field of external sources is known numerically. 

The present derivation is much Simpler since we assume the 
Earth to be represented by a medium of infinite conductivity 
from the center out to a radial distance b (see Figs Agate 
and a medium of zero conductivity from the radial distance b 
to the surface. 

In the third section we apply the results of this 
induction problem to obtain the magnetic field due to 
currents induced in the Earth by the arbitrary three- 
dimensional current system studied in Appendix Al.1. 

The fourth section summarizes the results obtained 
by Bostrém (1971) who used the image dipole method to 
obtain the field of currents induced by the three-dimensional 


current system discussed in Appendix Al.2. 
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A2.2 A simplification of Ashour's method 
ee OT nour $s method 


We assume the Earth to be composed of an infinitely 
conducting sphere of radius b Surrounded by a shell with 
zero conductivity and thickness a-b where a is the 
Lemeteeor thee larthc (sees rig. .A2,4¢)) awe antelaséume the 
external source to be a static current system external to 
the Earth such that the field can be calculated numerically 
everywhere inside and on the surface of the Earth. We can 
thus represent the scalar potential for the external field 


by a series of spherical harmonics as 


and the scalar potential for the field of the induced current 
as 
wales Cibei F' 
j = m gis imo 
V(r 9285209) = bo Dd Bay Py (eos @) (5% el MPR (im 
(A2.2) 
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n,m n=0 m=0 


The three components of the inducing and induced 


fields in the directions (r 285 29,) are obtained as follows: 
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=/ - F nv A pl (cos 8) () ei md (A223) 
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(ro18o, po) 







rg 
see 
Xx SURFACE OF 
SPHERE (o-= 00) 
Fig. A2.1 Diagram showing the spherical conductor at depth 


a-b and observation point (1 28559,) where the total 
field is calculated. The dashed line from the 
Origin to ro along the position vector . represents 
the path of integration for calculating the hori- 
zontal components of the induced field. 
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Since we assume the coefficients ner are known 
for the external field, the coefficients Et for the 
induced field can be determined by applying the proper 
boundary conditions at the surface of the infinitely con- 
ducting sphere. We have a choice of using one of the three 


following boundary conditions, 
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rae 
Pe) e M 
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or 0 ) 3740. OF Fp aa aie 30 
r_=b 


(A2.13) 


(see Stratton, 1941). 


Applying boundary condition A2.11 for the vertical fields 
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Hence, A2.6, A2.7 and A2.8 become 
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respectively. 

Expressions A2.15 through A2.17 for the components 
of the induced field are equivalent to those obtained by 
Ashour (1971) for the high frequency and infinite conducti- 
Perveiimiaceror a- thin spherical snell of radius 6b . The 
following derivation of the relationships between the inducing 
and induced fields at the surface of the Earth is essentially 
the same as that presented by Ashour (1971). 

Since we are interested in the value of the induced 


field at the surface of the Earth ar = a) , the components 
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of the induced field are represented in terms of components 
of the external field at r._=a. It can easily be shown 


from A2.3 and A2.15 that 


sp) (A2.18) 


where 


2 
0 Aa 


To obtain similar relationships between the horizontal 


components of Be and B' is made more difficult due to 
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be written a 


ee 1 = oe (A2.19) 
Bs and B can be represented by 
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where 
os Bar ,) m im 
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(A? .23) 
~1] r -(nt2) 
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The terms 4 and is, can be related to the external 
fields BS and BS in the same manner as for the vertical 


field By to give 
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Cr (4.05565) = (2) Bg (r.,0,4¢,) (A2. 26) 
and 
o,) = (2) BE (r.,0, 44.) (A2.27) 


The terms D, and D, can be represented by an integration 


of ee and ribs along the position vector of the 
observing coordinates (r 2852) from the origin to oe 


(see Fig. A2.1) as follows: 
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Therefore, since the external field can be evaluated 
numerically along this path of integration, the total in- 


duced field at Daten ee lie An be evaluated numerically using 
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Induction problem using Biot-Savart's law 
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(A2-31)) 


In the previous section equations A2.18, A2.30 


anid: wvd.3) 


surface of the Earth were derived. 


relating the induced and inducing fields on the 


In this section we 


substitute into these equations the inducing field equations 
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arbitrary three-dimensional current system in the magneto- 


Sphere and induced currents 
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We showed in appendix Al that for any external 
current’ Toop’ C’ {see Figie A2.1)° with tds representing an 
element of current along its path, the components of the 


fretd at (1 598525) in the coordinate directions (r 662%) 


are given by Al.21. For convenience, elements of dc are 


Oo 


Substituted in Al.21 to give three component equations 


e 7 r 
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We see from these equations that the vertical 
component and the first terms of the two horizontal compo- 


nents of the induced field can be obtained by substituting 


ror ts in Ag.32 through Az.3¢ and multiplying by the 
constant factor (2) with the appropriate sign. The main 
difficulty that arises in calculating the induced field is 
the evaluation of the integral term in both horizontal 
component equations A2.30 and A2.31. This uote Ost it Varies 
partially overcome by the fact that the integration over 
the source (line integral along C) is independent of the 


variable a thus enabling us to change the order of 


imiegration “(i.e for ‘the BS component) as follows, 
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Since the matrix A and differential line element ds are 
independent of ry» we can integrate directly over terms 


3 


such as Gait and re 7R° . Therefore, we have to find 


solutions to integrals of the form 
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The Avy element of matrix A_ represents the 
cosine of the angle i between position vectors rand ef 
and hence, can range in value from -1 to +1. Thus, the 
integrals—lI and I must be solved separately for the 


two cases, 
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resulting incthe following four integrals 
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The absolute value sign in A2.39 can be neglected since 


ee evey the interval of integration. The solutions 


of these integrals are listed below: 
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Case II Geos. 8. # 7)) 
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le = —— 
) a 


Since we must integrate terms containing I, and 


I over the entire current loop C (see equation A2.35), 


2 
the task of switching between case I and case II depending 
upon the value of cos X along the path of integration 
is very difficult due to the singularities in A2.40 and 
AZ.41 at cos } = +1. As a result of these singularities 


any numerical evaluation of the integral A2.35 is suspect 
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because of the change of solution for ] anger of from 
Case | to case Il ag cos) approaches +7] 

However, we can overcome this problem by noting 
tnat siege) ree nh 3 solutions given by A2.40 and A2.41 can 
be rewritten to remove the singularities at cos i = ap 


The solution for I, (A2.40) can be rewritten as follows: 
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It can be shown that A2.44 is valid for all values of 
COS A,eand, Tor .cos +4 =.+ sR reduces gopthesasolution 
Revtes? Mne Tirst term of the selutionefor I, (AZ 41h) 


can be rewritten in a similar manner, 
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‘ i cos? A - 1) ye KRA™ r) cos i 
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We must, however, consider two separate solutions in place 

of the second term because both numerator and denominator 
contain singularities. Since our main concern is to 

obtain proper solutions in the neighborhood of cos i = tas 
we choose two regions, namely cos X} < 0 and cos A> O 

for the determination of these separate solutions. For 


coS Xk < 0 the second term 
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remains unchanged. For cos A > 0 
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R + 3 - r cos A 
13 = log | ———] 
rid cos 1) 


a OhiG r{jod cos) | (A2.47) 
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a Pals A 
is cos i 


Therefore, 


a b- 
I, (cos A < 0) = I, + I, (A2.48) 
and 
_ a b+ 
I, feos A’ > 0) = I, + I, (A2.49) 


Equation A2.48 and A2.49 reduce to A2.43 for cos A = -]1 
and cos 4 = 1 respectively. We now have a set of solutions 
for integrals I, and I, valid for all values of cos A 


provided r_< ur Hence, A2.35 becomes 
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+ I,A35)ds, + (I,A, - rTyApp)d5] (A2.50a) 


15A,5)dS, - (1,A,3 + rTyAyp)455| (A2.50b) 


By substituting A2.50a and A2.50b into A2.30 and 
A2.31 respectively we obtain a complete solution for the 
induced field due to an arbitrary current loop evaluated 
by Biot-Savart's law. Using Al.20 and A2.18 we can represent 


this complete solution by the following equation: 
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(ey = -l, Eo .3 = +) (AD 51) 


where dc is obtained merely by substituting t for 


Be in the matrix dc and des. are the elements of the 


matrix 
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A2.4 Induction problem using the image dipole method 


Bostrom (1971) expanded equation Al.25 (see 
Appendix Al) in spherical harmonics in order to obtain a 
solution for the induced field associated with the current 
loops shown in Fig. Al.2 external to a spherical conductor 
mein iffite conauctivity (see Fig. Ace l). BY satistying 
the boundary condition for the vertical field (see eq. 
A2.11) he was able to show that the scalar magnetic poten- 
tial due to the induced currents on the surface of the 


spherical conductor of radius b is given by the equation 
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(r= b“/r) 


The expression for the magnetic field due to the 


induced currents is obtained by replacing dm, aie ou e 


‘ 
R by d°m, Chita and = in eq. Al.28. The evaluation 


of the partial derivatives is much more complicated than 


for the external field due to the added term in both R. 


and R, . TO remove some of the difficulties with evaluating 


these partial derivatives, Bostrom derived the following 


expressions: 
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(A2.53) 
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APPENDIX A3 THE METHOD OF GAUSSIAN INTEGRATION 


Since the integrals presented in Chapter II and 
Appendices Al and A2 cannot be integrated analytically we 
must resort to numerical methods for evaluating such inte- 
grals. The method of numerical integration we found to be 
the most suitable in terms of accuracy and efficiency is 
the Gaussian method. Gauss' formula for integration is 


given as follows: 


— 


n 
F(xdax boy) w, f(x.) +E (A3.1) 
(21 : : 
-] 
where Xs is the jth zero of the Legendre polynomial 
P(x) and 
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Numerous tables of X and Ws for various n- are 


readily available (see Abramowitz and Stegun (1965)). The 


error term EY given by 
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4 
2ntl 
For 4 peetegre ls th 662M) (yy (-Lsn<t) (A3.2) 


(2n+1)[ (2n)! ] 


can be approximated provided the anth 


derivative of f(x) 
can be evaluated. This is, however, not possible for the 
functions we must integrate. The determination of the 
error is made even more difficult due to the fact that we 
have to evaluate double and triple integrals. A discussion 
of E, for double integrals is given by Ralston (1965). 
Another factor affecting the error term EY is the change 
of variable necessary to make the interval of integration 
-I1 <x <1. For example if we have an integral with an 


interval of.integration- a <1 <b, a change of variable 


can be made to give 
] 
_  b-a b-a 
dew Eke = ie fla = aa (x+1) ]dx (A893 ) 
-] 


The effect of this change of variable on the error term is 


pie) die Letoneed 565k 


a multiplication factor [(b-a)/2] 
Thus, as long as the interval b-a is less than 2, the error 


will be reduced. If not it may increase significantly. 
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since the error term EY could not be evaluated 
for the integrals discussed in Chapter II and Appendices 
Al and A2, we subjected them to various degrees of Gaussian 
integration until convergence of the numerical results for 
each integral was obtained. As a further test for conver- 
gence we subdivided the intervals of integration, thus being 
assured of accurate numerical results. 

The final and most significant test was performed 
by integrating over model current systems using both methods 
(see Appendices Al and A2) of evaluating the associated 
magnetic fields. We found that the numerical results using 
these two methods were surprisingly similar for all model 
current systems integrated. Table A3.1 illustrates this 
Similarity for an E-W current system. The induced field 
is also included in these results. It is interesting to 
note that the method derived in this thesis (Biot-Savart 
method) requires much less computer time (a factor of 6 to 
10) than the magnetic dipole method as derived by Bonnevier 
et al. (1970) and Bostrom (1971). This fact is highly 
Significant when considering large current systems or 
comparison of model field calculations with a large number 


of magnetic observations. 
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APPENDIX A4 INSTRUMENTATION 


The basic system used to measure and record 
variations of the Earth's magnetic field at each field 
station is shown in Fig. A4.1. Each system included a 
three-component fluxgate magnetometer, analogue-to-digital 
converter, digital tape recorder, WWVB radio receiver and 
the necessary external power supplies and regulators. 
Since geomagnetic storms and related ionospheric distur- 
bances have little or no effect on the reception of WWVB 
time transmissions (60 kHz carrier frequency), the recording 
of the WWVB signal for two minutes at 7.5 hr intervals 
provided excellent timing for the magnetic records. 

All except one of the magnetometers used in this 
experiment were of the type designed by Dr. P. H. 

Serson of the Geomagnetic Division of the Department of 
Energy, Mines and Resources (see Fig. A4.2). The design 
specifications and results of tests performed on the 
magnetometer under normal field conditions have been 
presented in detail by Trigg et al. (1970). The Sharpe 
(MFO-3) fluxgate magnetometer used at Leduc has similar 
specifications except for better long term stability 

eos ly/24 hrs at constant temperature), a smaller temperature 


coefficient (+ 0.5y/°C) and a broader frequency response 
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A4.2 Block diagram of the fluxgate magnetometer. 
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(DyG. to 55 Hz, <Sdb).. The theoretical frequency response 
for the magnetometer is shown jn Pugs “Ads .o, Lt Can 

be seen from this figure that the phase and gain curves 
are essentially flat with the phase shift being zero from 
Oeeta ) 05025 He , the frequency interval used in this 
thesis. The sensitivity, resolution and dynamic range of 
the magnetometer are 1 volt/100y , ly and +1000y (from 
baseline) respectively. The Outputs Wrrom the X , ¥ and 
Z amplifiers are linear over most of the dynamic range 
(+10v) (see Fig. A4.4). 

The stability and drift tests completed by Trigg 
et al. (1970) showed that fluctuations of the Y baseline 
(up to 20y) were much larger than expected and could not be 
explained by temperature drifts. Also a temperature 
coefficient of 10y/°C was necessary to explain the 
variations of the Z baseline. The value of the temper- 
ature coefficient given by the manufacturer of the magneto- 
WetSers 1S. < 1y/°C 

The sensing head (Fig. A4.3) was mounted on an 
aluminum pole set in concrete at the bottom of a plastic 
container (waste receptable). To reduce diurnal tempera- 
ture variations the container was buried in the ground, the 


top being ~ 6 to 10 in. below the surface. Although 
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Theoretical frequency response curves for the 
magnetometer (after Trigg et al. (1970)). 
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A4.4 Gain curve of the magnetometer for the lower 
portion of the dynamic range. 
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fluctuations in temperature were not determined, Ievoog et. al. 
(1970) suggested that burying the sensing head in the ground 
can reduce the daily temperature variations to one or two 
degrees centigrade. 

The Redcor 720 A-D converters used in this 
experiment were modified to operate at various stepping 
rates, thereby allowing magnetic data and WWVB time signals 
to be recorded with the same digital system. After each 
series of 7 blocks of data (~ 65.5 min/block) the output 
of magnetic data from the multiplexer was interrupted for 
123 sec while WWVB time signal was being recorded (Fig. 
A4.5). To permit accurate recording of WWVB the stepping 
rate was changed from 1.56/sec (sample interval of 1.92 
sec/component) to 50.0/sec. These stepping rates were 
controlled by a 400 Hz crystal oscillator (with a maximum 
deviation from nominal frequency of + 0.005 per cent) and 
an associated divider circuit. Hence the maximum variation 
in the nominal 7 hr 38 min interval between successive 
WWVB recordings was %~ 2 sec. 

The cyclic sampling order (HDZHDZ...) of equal 
sample intervals in the multiplexer was controlled by the 
sequencer. The 12 bit A-D converter (the least significant 


bit was used as a flag to identify errors in the Sampling 
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order) output was formatted into two byte samples (each 

byte being 6 bits long) and recorded on a seven-track 
Peripheral Equipment incremental write tape deck. Each 
block contained 6144 samples and represented 65 min 32.2 sec 
of magnetic data or 122.8 sec of WWVB time signal. Since 
the dynamic range was + 1000y , the least significant bit 
represented 0.976y (2000y/2!!) which is approximately the 
Same as the resolution of the magnetometers (+ ly). 

The component outputs from the magnetometer were 
passed through two stage aliasing filters to remove frequency 
components above the Nyquist frequency (0.26 Hz) (see Fig. 
A4.6). The filters also served as buffers to isolate the 
multiplexer from the magnetometer by limiting the voltage 
input to less than 12 volts (+). 

Because of the loss of significant amounts of 
data during the 1969 field season due to power failures, 
an automatic restart module was designed and installed at 
each station. Immediately following any power failures 
the restart process for the digital system was automatically 
carried out. Also the WWVB time signal was recorded 
immediately following the power failure to increase the 


timing reliability. 
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APPENDIX AS INTERPRETATION OF MODEL CURRENT SISTEMS 
WITH A HYPOTHETICAL LINE OF STATIONS 


In this appendix we show examples of reproduction 
of theoretical profiles obtained by assuming that model 
current systems with different latitudinal current distribu- 
tions were observed by a hypothetical line of stations with 
the same distribution in latitude as our line of stations. 
Latitude plots were made using the values obtained from three 
different current systems, each positioned at different 
latitudes. Latitude profiles were then interpreted from 
these plots in the same manner as the actual data from our 
line of stations, thereby, allowing a comparison to be made 
with the original theoretical profiles. For this exercise 
we assumed all stations to be in operation. 

The Tirst*\ example is shown.in Fig. A5.1.\\Profile A 
is the theoretical profile whereas profiles B, C and D are 
the interpreted profiles for the model current system located 
at tie 70.0° and 72.5° respectively. The interpreted 
profiles for the model current system at these three different 
latitudes agree reasonably well with the theoretical profile. 
The value of H at 70° in profile C was inferred from the value 


Qf Heand 2 at the Southern anid northern’ stations. For the 
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current system at 72.5° the Interpretation was more difficult 
due to the lack of stations between 71.4° and 77.0°. In the 
present study the magnetic data from Baker Lake has proved 
to be invaluable for the interpretation of high latitude 
current systems even though the station is 18° east of our 
line. Again for profile D the interpretation was based on 
the southern and northern stations. 

PYGaethOse Snows the results foreaawide current 
System with an enhanced northern border. For this example 
the reproduction is very good for the current system being 
located at all latitudes shown. This is not however the 
case for the double current system shown in Fig. A5.3. 
Although profile B is in good agreement with the theoretical 
profile, the interpretation was more difficult when the 
current system was located at higher latitudes. In fact 
profile D completely misrepresents the actual current system. 
Thus, when the interpretation of a current system is dependent 
upon stations at lower or higher latitudes, such an inter- 
pretation may fail to show the more complex distributions 
in the current flow along the auroral oval (i.e. double 


current systems). 
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PYeEN LnouGn the resolution of our line of stations 
is difficult to ascertain since the interpretation between 
stations for one component may be accomplished with the aid 
of other components either or both north and south of this 
particular reqime, we feel the results of the foregoing 
exercise enhances our confidence in the interpretation of 


actual current systems presented in this thesis. 
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APPENDIX A6 THREE-DIMENSIONAL DIAGRAMS 


In this appendix we present a selection of three- 
dimensional (3-D) diagrams of some of the polar magnetic 
substorms discussed in Chapters IV and V. It should be 
noted that the sign of the H component has been reversed. 
Hence a westward electrojet has the appearance of a moun- 
tain or mountain range while an eastward electrojet has the 
appearance of a valley. The amplitude scales on these 
plots are different for all components. Therefore a 
comparison between the relative amplitudes of different 


components can not be made with these particular diagrams. 
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APPENDIX A7 MORE EXAMPLES OF MAGNETIC PERTURBATIONS 
PooVC VATED SWITH POLAR ELECTROJETS 


Several examples of polar magnetic substorm 
events and magnetic perturbations associated with east- 
ward electrojets are presented in this appendix. Two 
latitude profiles for each event are also included. The 
Kp indices for the days represented in this appendix are 


7¥sted tW Table A7.T% 
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Table A7.1 Kp indices for the days represented in 


this appendix 


Three-hour range indices 
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Fig. A7.1 Component magnetograms for June 19, 1970 
(0430-0800 UT). The existence of a substan- 
alee Ue curiat ion. at floMeand PF lCH and tie 
+Z perturbation at CAMB prior to 0445 UT are 
indicative of a westward electrojet between 
CAMB and FTSM. Note the large +D perturbation 
geri ates. 0600 Uda shis was. propab ly 
associated with a westward travelling surge. 


Fig. A7.2 Component magnetogram for June 19, 1970 
(0830-1200 UT). This is an example of a 
large bay type event commonly seen in the 
post midnight sector. 


Fig. A7.3 Latitude profiles for the events shown in 
Pigeerd st and Aree. eenoLce the simitarity 
between the H and D profiles at 0610 UT. 

The electrojet was probably orientated in 
the NW-SE direction at this time. We have 
observed this effect several times, especially 


in the evening sector. 


Fig. A7.4 Component magnetograms for June 30, 1970 
(0500-1000 UT). Note the buildup of an 
eastward electrojet prior to the substorm. 
This event resembles the event (Day 195) 
discussed in Chapter IV. 
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Component magnetograms for September 18, 1970 
(0600-1000 UT). Note the transition of the 

D component from - to + at all latitudes. 
This particular event was associated with 
Slowly drifting N-S aligned auroral arcs. 
This D component signature was observed on 
other occasions. 


Latitude profiles for the events shown in 
Fig. A7.4 and A7.5. Note the development 
of a new current system at the northern 
border of the pre-existing electrojet in 
each case. 


Component magnetograms for July 2, 1970 
(0615-0915 UT). This particular event was 
discussed in detail by Kisabeth and Rostoker 
(1971). 


Component magnetogram for July 2, 1970 
(0930-1300 UT). 


Latitude profiles for the events shown in 
Aidwuerwet and Az 8. -Protlies. Avani 6 
illustrate the decay of the northern border 
leaving the stable portion .of-the, electrojet 
to the south. Profiles C and D show the 
enhancement of the northern border of an 
electrojet in the post-midnight sector. 
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Component magnetograms for July 15, 1970 
(0830-1200 UT). Note the sharp onset and 
the associated +H perturbations at lower 
latitudes. 


Component magnetograms for July 12, 1970 
(0600-1100 UT). This is an example of a 
bay type event of long duration. The long 
period perturbations superimposed upon 
this bay represent variations in the 
westward electrojet intensity. 


Latitude profiles for the events shown in 
Fig. A7.10 and A7.11. Note the northward 
motion of the Z crossover in profiles A 
and: &.- =Profites-€ and- D-depict the 
enhancement of the northern border. 


Component magnetograms for September 19, 1970 
(0440-1140 UT). Shown here are examples of 
substorms occurring in both the evening and 
midnight sectors. Note the long magnetic 

bay event (». 0630-1130 UT) with a substorm 
Superimposed upon it. 


Latitude profiles for the events shown in 
Fig. A7.13. Note the development of a double 
current system in the evening sector as well 
as in the midnight sector. 
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Component magnetograms for July 8-9, 1970 
(2300-0230 UT). This is an example of an 
SSC recorded with our line of stations. 


Component magnetograms for July 9, 1970 
(0230-0630 UT). Note the buildup of a 
large eastward electrojet along with a 
later development of a westward electrojet 
to the north of the eastward electrojet. 
Also note the intense low latitude (x 62°) 
magnetic substorm onsetting at » 0545 UT. 
ineoh index TOY tie interval of. time 
including the onset of this substorm is 
7-, thus indicating an expanded auroral 
Oval. 


Latitude profiles for the events shown in 
Fig. A7.15 and A7.16. Note the development 
of an eastward and westward electrojet 
following the SSC as illustrated in profiles 
A and B. Profiles C and D show the develop- 
ment of a large westward electrojet north 

of a pre-existing electrojet. 
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